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ABSTRACT 

VLA 1.4 GHz (cr - 0.012 mJy) and MIPS 24 and 70 ^m (a ~ 0.02 and 1.7mJy, respectively) 
observations covering the 2 deg 2 COSMOS field are combined with an extensive multi-wavelength 
data set to study the evolution of the IR-radio relation at intermediate and high redshift. With 
~4500 sources - of which ^30% have spectroscopic redshifts - the current sample is significantly 
larger than previous ones used for the same purpose. Both monochromatic IR/radio flux ratios (^24 
& g7o ), as well as the ratio of the total IR and the 1.4 GHz luminosity (<?tir) are used as indicators 
for the IR/radio properties of star forming galaxies and AGN. 

Using a sample jointly selected at IR and radio wavelengths in order to reduce selection biases, we 
provide firm support for previous findings that the IR-radio relation remains unchanged out to at 
least z ~ 1.4. Moreover, based on data from ~150 objects we also find that the local relation likely 
still holds at z 6 [2.5, 5]. At redshift z < 1.4 we observe that radio-quiet AGN populate the locus 
of the IR-radio relation in similar numbers as star forming sources. In our analysis we employ the 
methods of survival analysis in order to ensure a statistically sound treatment of flux limits arising 
from non-detections. We determine the observed shift in average IR/radio properties of IR- and radio- 
selected populations and show that it can reconcile apparently discrepant measurements presented in 
the literature. Finally, we also investigate variations of the IR/radio ratio with IR and radio luminosity 
and find that it hardly varies with IR luminosity but is a decreasing function of radio luminosity. 
Subject headings: cosmology: observations - galaxies: active - galaxies: evolution - galaxies: high- 
redshift - infrared: galaxies - radio continuum: galaxies - surveys 
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1. INTRODUCTION 

The global infrared (IR) and radio emission are 
tightly and virtually linearly correlated in a broad va- 
riety of star forming (SF) systems (see iHeloul 119911 : 
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iCondonl 119921 : lYun et all 120011 and references therein), 
thus defining what is known as the ' IR-radio rela- 
tion'. Studies in the nearby universe have shown that 
not only late-type galaxie s (Dickcv & Salp eterl 119841: 



Helou et alJll985HWunderlich et alJH987HHumrnel et alJ 
19881) ranging from normal spirals to the most vigorously 



star forming ultra-luminous infrared galaxies (UL IRGs; 
e.g. iSanders fe Mirabell [l99l iBressan et all [20021) fol- 
low the relation, but also many early-type galaxies 
with low- level star form ation (|Wrobel fc Heesch"enlll988t 
iBally fc Thronsonl 119891 ) . as well as int eracting systems 
of m ixed morphological composition (|Domingue et all 
120051) . 

Following first indications of the correlat ion in ground- 
based observations at lO^ m and 1.4 GHz (jvan der Kruitl 
119731 : [Condon et al~lll982l ) the ubiquity and tightness of 
the IR-radio relation became fully appreciated during 
the analysis of the combination of data from the Very 
Large Array (VLA) and the Infrared Astronomical Satel- 
lite (IRAS) which measured the far-I R (FIR) properties 
of ~20,000 galaxies at z 5 , 0-15 (e.g. iDickev fc Sahoetei 
19841 Ide Jong et alJ [1985Y iHelou et alJ 11985b lYun et al " 
20011 ). IR observations of sources at higher red- 



shifts became available with the advent of the Infrared 
Space Observatory (ISO) and, in recent years, with the 
Spitzer Space Telecope. They have provided increas- 
ing evidence that t he locally observed correlation likely 
holds until z ~ 1 (|Garrettl l200l lAppleton etafl 120041 : 
IFraver et alJ 120061 ) and that the linearity of the cor- 
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relation is maintained as far back as z ~ 3 although 
the slope 17 may change, especially for sub-mm galax- 



ies dKovacs et al.ll2006HVlahakis et a Tll2l)07t lSaiina et al 



20081: IMurphv et alJ l2009bL but see also iBeelen et al 
20061: llbar et al.1 120081)" The statistical significance of 



these high-z studies, however, is still low because so far 
the number of sources detected at z > 0.5 is limited. 

The very tightness of the IR-radio relation (the intrin- 
sic dispersion in the local galaxy population is less than a 
factor 1.5), combined with the fact that it spans five or- 
ders of magnitude in bolometric luminosity has provided 
a useful tool for numerous astrophysical ap plications and 
motiv ated continued study of the relation. Brcssa n'et al.1 
(2002), for example, use the measured IR/radio flux ra- 
tio to determine the evolutionary stage of systems un- 
dergoing a starburst, depending on whether they show 
excess IR or radio emission with respect to the aver- 
age locus of star forming galaxies (SFGs). The con- 
sideration of radio-excess outliers to the IR-radio rela- 
tion also has been used to select galaxies in which an 
active galactic nucleus (AGN) rather than star forma- 
tion is the dominant source of IR and radio emission 
(e.g. IDonlev et ai]|2005t IPark et al.ll200l . At high red- 
shift the IR-radio relation has been used to compute dis- 
tance est imates of sub-mm galaxies lacking optic al coun- 
terparts (|Carilli fc Yunlll999l:lDunne et aljboOOj) . as well 
as to estimate the contribution of SFGs to t he cosmic 
radio background (Haarsm a fc Part ridge 1998|). 
In the absence of AGN activity both (F)IR and radio flux 
measurements are in principle unbiased tracers of star 
formation since the radiation in these spectral regions 
is not attenuated by dust as opposed to the often heav- 
ily obscured emission at ultraviolet (UV) or optical to 
near-IR (NIR) wavelengths. The observation that a cor- 
relation between the IR and radio exists implies that the 
thermal IR dust emission and the radio continuum flux 
which is a frequency-dependent mixture of thermal emis- 
sion and a non-thermal synchrotron component - have a 
common origin. Soon after the discovery of the IR-radio 
relation the birth and demise of massi ve (>5M fl ) stars 
was i dentified as it s likely cause (see lHarwit fc Pacinil 
119751 : [Condonlll992l . and references therein). Deviations 
from exact proportionality of the IR and radio emission 
as well as changes in the IR-radio relation with cosmic 
time thus potentially imply changes in the mechanisms 
steering star formation, or that IR and radio emission are 
not equally good proxies of star formation under all cir- 
cumstances. However, several decades after the discovery 
of the relation many details concerning the physical pro- 
cesses which shape it still need to be settled. 
A number of models to explain the global correlation 
between the integrated IR and radio fluxes of SFGs have 
been advanced: they include the calorimeter model (|Volk 



1989 ] IPohl[l99l ILisenfeld et alJlToM iThompson etal 

2004 I200I 



„ , the 'optically thin' scenario (|Helou fc B icav 
19931) and the 'linearity by conspiracy' picture of iBell 
(2003). Models focusing on the magnetohydrodynamic 
state of the interstellar medium (ISM) are often con- 
structed in order to explain the correlation on both 
kpc scales within galaxies as we l l as globally (e.g. 



iBettens et all Il993t iNiklas fc Beck! Il997t iGroves et al.1 

17 By 'slope' we mean the slope of the correlation in a plot of 
flux vs. flux with linearly scaled axes. 



120031: iMurgia et all 120051 ). Many of these fare well in 
reproducing multiple aspects of the correlation but ei- 
ther cannot satisfy all observational constraints or make 
predictions which still await confirmation. 
On the observational side further insight into the (astro- 
physical processes shaping the IR-radio relation has 
been gained by studying the IR-radio rel ation on small 
scale s in resolved nearby galaxies (e.g. iBeck fc Gollal 



[19881: iBicav fc Heloul [19901: IMurphv et all 120081 Dumas 
et al., in prep.) which revealed that IR/radio proper- 
ties vary inside spiral galaxies or by investigating the 
impact of environmental effects (I Miller fc Owenl 120011 : 
iReddv fc Yunl 12004 IMurphv et alll2009al ). Another im- 
portant finding was that the linearity of the relation is 
maintaine d with both the therm al and non-thermal radio 
emission (jPrice fc Durid fl992l ) taken separ ately, as well 
as for the warm and cold dust components (jPierini et al.1 
l2003h . 

As a q uantitative measure of the correlation iHelou et al.1 
( 1985) introduced the logarithmic flux density ratio 



<?fir = log 



FIR 



3.75 x 10 12 Wm" 2 



/ S V (1A GHz) \ 
' ° g ^Wm^Hz- 1 J 



(1) 

Here -Ffir is the rest frame FIR flux - traditionally com- 
puted as a linear combination of measurements in the 
IRAS 60 and 100 /im bands under the assumption of a 
typical dust temperature of ~30K - and 3.75xl0 12 Hz 
is the 'central' frequency of the FIR-window (42.5- 
122.5 /im). S V (1A GHz) is the monochromatic rest frame 
1.4 GHz flux density. 

Strong improvements of the sensitivity of IR observa- 
tories over the last decade, in particular in the mid- 
IR (MIR) with the MIPS photometer on Spitzer, have 
promoted the use of monochromatic flux ratios, e.g. at 
24 or 70 urn (e ,g. lAppleton et al.l l200l llbar e tal] 120081: 
iSevmour et al.ll2009l) : " 



, /5„(24[70]/im)\ 

924 [70]= log ^ (14GHz) ) , 



(2) 



where both flux densities are specified in units of 
Wm~ 2 Hz -1 . This approach is convenient for evolution- 
ary studies of the IR-radio relation at intermediate and 
high rcdshift as it avoids the computation of IR luminosi- 
ties from a single IR flux measurement, usually made at 
24 fim due to the high sensitivity of the according MIPS 
filter. However, it still requires the computation of K- 
corrections, if the flux densities in equation ^ are to be 
given in the rest frame. 

Since current radio surveys at centimeter wavelengths 
are generally significantly shallower than the MIR pho- 
tometry several authors have carried out radio stack- 
ing experiments (iBovle et al.1 120071 IBeswick et~al]|2008t 
IGarn fc Alexanderl I2009T with the aim of determining 
the IR/radio properties of the faint IR population. This 
is of particular interest in view of the potential for de- 
tecting changes in the IR-radio relation brought about 
by relativistic cooling of cosmic ray electrons by in- 
verse Compton scattering off photons from the cosmic 
microwave background (CMB). This cooling may over- 
whelm - at least at low decimeter frequencies - the syn- 
chrotron losses in t he ISM of normal galaxies at z > 0.5 
(|Carilli et al.l 12008) if the energy density of their mag- 
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netic fields is similar in strength to that ty pically mea - 
sured in spiral arms (Ub ~ 10 12 ergcm~ - 3 ; Beck' 2005). 
Alas, the stacking results have produced strongly dis- 
crepant results, both among different cosmological sur- 
vey fields and with respect to recent studies of sources 
directly detected in both the IR and radio. An additional 
concern of particular relevance in evolutionary studies is 
the bias that is introduced by constructing samples from 
different selection bands in flux limited surveys. 
In this paper we focus on the impact of selection biases 
and apply the statistical technique of survival analysis to 
our data which permits the inclusion of constraints from 
flux limits in the study. This approach is a significant 
improvement over limiting the sample to those sources 
which are detected at both IR and radio wavelengths. 
These topics are discussed in §[S] Our data sets contain 
an unprecedented amount and quality o f information 
gathe red as part of the COSMOS survey (jScoville et al.l 
2007). This is reflected by a sample which contains 
significantly more sources in the redshift range z > 0.5, 
in which the sources in previous studies began to taper 
out. Moreover, a large fraction (33%) of spectroscopi- 
cally determined redshifts, which are supplemented by 
accurate photometric redshift estimates, in combination 
with flux constraints at both 24 and 70 /im for each 
source in the sa mple make for impro ved estimates of 
IR luminosities (jMurphv et al.l I2009U e.g.; Kartaltepe 
et al.j subm.). They are an important prerequisite for 
placing accurate constraints on the evolution of the 
IR-radio at high redshift. We introduce our data and 
the subsequently analyzed samples in §[5] For those 
interested, a detailed review of the different catalogs 
and images used in the paper and the band-merging 
of this information into a rich multi-wavelength data 
is provided in Appendices [A] and [B] §§[3] and H] deal 
with methodological considerations applying to the 
identification of star forming sources and AGN galaxies 
as well as to the derivation of IR luminosities. Following 
the section on biasing and our treatment of flux limits 
in §[5j the bulk of our analysis is then presented in §[6l 
Our most important results are: 

(a) the constancy of the IR-radio relation as 
parametrized by the flux ratios (724, Q70 and <7tir 
out to z ~ 1.4, as well as for a sub-sample of high-z 
sources at 2.5 < z < 5, 

(b) the identification of selection biases as a potential 
explanation for discrepant average IR/radio flux ratios 
measured in previous studies, and 

(c) the observation that over the last 10 billion years 
the distribution of IR/radio ratios of optically selected, 
radio-quiet 18 AGN has been very similar to that of star 
forming galaxies 

We discuss and summarize these findings in §§0 and [51 
respectively. 

Throughout this article we adopt the WMAP-5 cos- 
mology defined by £l m — 0.258, il\ + & m = 1 and 
a present-day Hubbl e parameter of 71.9 kms _1 Mpc _1 
(iDunklev et al.ll2009fl Magnitudes are given in the AB- 
svstem of lOkd (fl974T ) unless the opposite is explicitly 

18 Objects referred to as 'radio-quiet AGN' in the rest of the 
paper are understood to be sources in which the radio emission 
from the AGN does not contribute significantly to the total energy 
emitted at 1.4 GHz. 



stated and we henceforth drop the subscript 'AB' in the 
text. 

2. DATA SETS AND SAMPLES 

In order to study biases arising from the selection of 
sources at either IR or radio wavelengths we have con- 
structed both a radio- and an IR-selected sample of COS- 
MOS galaxies. All observations and associated data sets, 
as well as the band-merging procedures used to identify 
counterparts from radio to X-r ay wavel engths are de- 
scribed in detail in Appendices IA.lflA.5l and Appendix 
[Bl respectively. Here we only briefly present the primary 
data sets fij l2.ll) and review the most important proper- 
ties of our radio-selected and IR-selected samples (§ !2.3p . 
This section also contains a summary of the ancillary 
multi-wavelength photometry and the redshift informa- 
tion which is available for our sources (§[22} • 

2.1. Radio and IR Data 

Th e VLA-COSMOS Project (ISchinnerer et al.l l200i 
120071 ) has imaged the COSMOS field at 1.4 GHz with the 
VLA to a mean sensitivity of <0.01 (~0.04) mJy/beam 
at the centre (edge) of the field. As the basis of our subse- 
quent analysis we use the VLA-COSMOS 'Joint' Catalog 
(Schinnerer et al. 2009, subm.) which contains ^2900 
sources detected with S/N > 5. Nearly 50% of these are 
resolved at a resolution (FWHM of synthesized beam) 
of 2.5" x 2.5". Flux measurements were carried out with 
AIPS (Astronomical Image Processing System; iGreisenl 
120031) and have been corrected for bandwidth smearing 
in the case of the unresolved radio sources. Their errors 
are generally ~17% of the flux value. 
IR data at 24 and 70 imp, we re taken by the S-COSMOS 
Survey ([Sanders et al.ll2007| ) using MIPS on Spitzer. 
At 24 /im the FWHM of the PSF is 5.8" and the aver- 
age la sensitivity ^0.018 mJy over a l arge fraction of 
the i maged area. Flux measurements (|LeFloc'h et al.l 
20091) in the dee P and crowded MIPS 24 /im image were 
performe d with the PSF-fitting algorithm DAOPHOT 
( Ste tson] [19871 ) which can simultaneously fit and hence 
de-blend multiple sources. In the following we use a flux 
limited catalog of COSMOS 24 /im sources which is re- 
stricted to the range 5V (24 fim) > 0.06 mJy. Typical flux 
uncertainties at 24 /jm are ~8%. 

MIPS 70 /im observations of the COSMOS field were 
carried out in parallel with the 24 /im imaging. They 
have an average la point source noise of 1.7 mJy and 
a resolution of 18.6". Our 70 /im source list includes 
detections down to S/N = 3 and was compiled as de- 
scribed in iFraver et al.l (|2009f ) . Fluxes were derived us- 
ing the APEX (Astronom ical Point source Extraction 
iMakovoz fc Mar lcau 2005) peak fitting algorithm; their 
average uncertainty is ~16%. 

In Fig. Q] we show the minimum total IR (TIR; 8- 
1000 /im) luminosity that is detectable as a function of 
redshift given the sensitivity of the 1.4 GHz (converted 
to an IR measurement assuming the local IR-radio re- 
lation) and 24/70 /im data. Different IR SED template 
libraries (see colour scheme in the lower right corner) lead 
to somewhat different predictions but it is clear that the 
1.4 GHz and 70 /im surveys have matching depths while 
the 24 /im observations are about seven times deeper. A 
similar sampling of the IR luminosity function is achieved 
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in all three bands if a 24 /tin flux limit of approx. 0.3 mJy 
is assumed. We therefore limit our 24 /tm catalog to the 
range S J/ (24/xm) > 0.3 mJy when we construct our IR- 
selected sample but allow fainter counterparts of 1.4 GHz 
sources to be included in the radio-selected sample. The 
use of these flux limited samples has the immediate con- 
sequence that we only detect the brightest ULIRGs at 
z > 1.5 while the average luminosity of our sources is 
much lower at, e.g. z ~ 0.5, where most sources belong 
to the LIRG class. 

2.2. Ancillary COSMOS Data 

Optical data and photometric red shifts are taken from 
the COSMOS photometry catalog of lllbert et all l|2009al) 
which lists more than 600,000 COSMOS galaxies with 
i + < 26 detected in a region roughly contiguous with the 
area covered by the VLA-COSMOS survey. The wave- 
length range covered by these observations (30 broad, 
medium and narrow band filters) extends all the way 
from the UV at 1550 A to the MIR at 8 fj,m. iCapak et all 
(|2007l I2008D provide a complete description of these ob- 
servations. 

Spectroscopic data has been gathered for more than 
20,000 sources in the COSMOS field, e.g. by the zCOS- 
MOS survey (ILillv et all l2007h and SDSS (lYork et al.l 
120001 ), or in Magellan/IMACS and Keck/Deimos follow- 
up o bservations dedicated to specific (clas ses of) sources 
(e.g. lTrump et alJ2007t [Trump et alll2009t Kartaltepe et 
al., in prep.; Salvato et al., in prep). If a reliable spectro- 
scopic redshift is available it is favoured over the photo- 
metric redshift estimate. The choice of the best possible 
distance measurement for our radio and IR sources is 
described in detail in 5 IA.5I 

The .XMM-Newton COSMOS Su rvey (IHasinger et al.l 
120071: ICappelluti et al.l 120071 120091) has detected a total 
of 1887 bright (>2xl0" 15 erg cm" 2 s" 1 in the 0.5-10 keV 
band) X-ray sources over 90% (1.92 deg 2 ) of the COS- 
MOS field. A large fraction of these are associated with 
AGN and hence provide a means of identifying AGN- 
powered radio and IR sources in our sample which is 
complementary to our primary classification scheme in- 
troduced in §[3] For our subsequent analysis w e rely 
on a list of XMM sources (|Cappelluti et all 120091 ) with 
unique and secure optical counterparts (see Brusa et al., 
in prep.) an d SEP fits to t he UV to MIR photometry 
performed bv lSalvato et al.l (|2009D . 

2.3. Description of the Samples 

Due to the differing characteristics (resolution, astro- 
metric accuracy) of the radio and IR data the determina- 
tion of counterparts at other wavelengths differed some- 
what for the radio- and the IR-selected sample. Which 
candidate counterparts are incorporated in the final sam- 
ple and which are rejected is determined by the goals 
of this study: in our case it is more important to se- 
lect objects with the cleanest possible radio and IR flux 
measurements rather than having a statistically complete 
sample. The details of the band-merging between the IR 
and the radio catalogs and the subsequent exclusion of 
ambiguous counterparts are presented in detail in Ap- 
pendix [B] Here we summarize the most important prop- 
erties of the radio- and IR-selected samples. 



2.3.1. The Radio-Selected Sample 

Based on the positions of >5er 1.4 GHz detections in 
the VLA-COSMOS Deep Project image we searched for 
IR counterparts in the S-COSMOS 24 and 70 fim cata- 
logs which have S/N > 3. Counterparts were determined 
by direct positional matching of radio and IR coordinates 
with search radii corresponding to approx. FWHM/3 of 
the IR PSFs of the respective MIPS bands. If no coun- 
terpart was found, a 3 a point-source detection limit was 
determined based on the corresponding uncertainty im- 
ages. Radio sources with ambiguous IR counterparts - 
i.e. in the presence of more than one potential counter- 
part or if the counterpart had not been uniquely assigned 
to a single radio source - have been excluded from the 
analysis of the paper. The match with the COSMOS 
multi- wavelength and spectroscopy catalogs provides dis- 
tance estimates for 73% of the radio-selected sample as 
well as photometry from the UV to the MIR which is 
used to separate galaxies dominated by star formation 
or AGN emission (see §[3]). In the upper panel of Fig. 
[2] we show histograms of the separation between radio 
source positions and the location of the optical and IR 
counterparts. Note that the distance is normalized by 
the width of the broader PSF of the two involved bands. 
Fig. [3] shows the 24 and 70 /tm flux distribution of the 
radio sources, including information on whether the flux 
constraint is a well-defined measurement or an upper flux 
limit. 

Fig. [4] (left-hand panel) explicitly shows how the fraction 
of sources that have a directly detected counterpart in ei- 
ther or both of the MIPS bands or only upper flux limits 
changes as a function of redshift (38% of the redshifts 
are spectroscopically, 64% photometrically determined). 

2.3.2. The IR-Selected Sample 

The IR-selected sample is based on sources listed in 
the S-COSMOS 24 /im catalog that satisfy the criterion 
S'„(24/Ltm) > 0.3 mJy. This criterion ensures that the IR- 
selected sample is well matched to the 70 fim and 1.4 GHz 
data as far as the sampling of the IR luminosity function 
is concerned. To reduce the likelihood of false identifi- 
cations due to the positional uncertainty of the 24 /mi 
sources we searched for IRAC counterparts, the posi- 
tions of which were used as a prior in the subsequent 
band-merging with the other wavelengths. If no IRAC 
counterpart was available we also admitted unambigu- 
ous matches with optical sources. 70 /jm counterparts to 
the 24 /im sources with S/N > 3 were determined and 
validated following exactly the same approach as in the 
radio-selected sample. For those 24 /im sources which 
did not already have a known radio counterpart (deter- 
mined in the construction of the radio-selected sample) 
we checked whether they are associated with a counter- 
part having S/N > 3. All new detections satisfying this 
criterion were then added to the list of radio counter- 
parts with S/N > 5 that were already known from the 
construction of the radio-selected sample. 24 /urn sources 
that are undetected at 70 /jm and/or 1.4 GHz are as- 
signed 3 a upper flux bounds. 

The distributions of separations between all 24 /im 
sources and their counterparts in the optical, far-IR and 
1.4 GHz maps are given in Fig. [2jb. The contribution 
of flux limits and well-defined flux measurements as a 
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function of flux and S/N at 70 fim and 1.4 GHz is illus- 
trated in Fig. [5] Finally, in the right-hand panel of Fig. 
[4] we show at different redshifts which fraction of the 
IR-selected sample has direct detections or upper flux 
density limits at 70/im and/or 1.4 GHz. Spectroscopic or 
photometric redshift measurements are available for 80% 
of the objects in the IR-selected sample. The remaining 
sources are either not bright enough for spectroscopy or 
have flux information in too few bands to derive a pho- 
tometric redshift based on SED fitting. 

2.3.3. The Jointly Selected Sample 

The jointly selected sample is the union of the radio- 
and IR-selected samples presented in §§ 12.3.11 and 12.3.21 
As such it contains 6863 sources: 1560 sources that are 
only detected at 1.4 GHz, 3960 sources that are only de- 
tected at 24 fim and, finally, 1341 sources which are se- 
lected at both wavelengths. In Table[T]we summarize the 
available redshift information for the jointly radio- and 
IR-selected sample, as well as separately for the radio- 
and IR-selected samples. 

In Fig. [6] the IR and radio fluxes of our sources are 
compared. The colour coding of the data points distin- 
guishes three kinds of sources: in black those which have 
entered both the 1.4 GHz catalog as well as the 24 fim 
catalog (restricted to sources with flux density larger 
than 0.3mJy), in green 1.4GHz sources without coun- 
terpart in the 24 /im catalog and in orange those 24 fim- 
detected sources which do not have a counterpart in the 
VLA-COSMOS Joint catalog. The sources from these 
three different categories have been added to the plot 
in random order to prevent that the symbols of the ini- 
tially plotted category are systematically hidden by the 
successively overplotted data in common regions of flux 
space. Fig. [6]c confronts the fluxes in the two selec- 
tion bands; the empty rectangle in the lower left corner 
of this panel reflects the selection criteria at 1.4 GHz and 
24 /im. Since the 24 /im catalog is flux limited, essentially 
all upper 24 fim flux limits lie at or below the critical 
flux threshold; upper 1.4 GHz flux limits for undetected 
24 fim sources on the other hand are also encountered 
at higher 1.4 GHz flux values than the sharp cut-off at 
5^(1.4 GHz) « 0.05 mJy because the radio-catalog was 
constructed using a S/N criterion. Note that the region 
where both the 1.4 GHz and 24 fim flux density clearly 
exceed the respective selection thresholds contains some 
sources which are not included in both the catalog of 
24 fim and that of 1.4 GHz detections (cf. orange and 
green symbols in the area where 5„(1.4 GHz) > 0.1 mJy 
and S'„(24 fim) > 0.3 mJy). Two reasons can be responsi- 
ble for this: (a) minor incompleteness of the catalogs, or 
(b ) spatial variations in the background noise which, at a 
given flux, lead to certain sources not being detected at 
the significance level required for inclusion in the original 
source list. 

3. IDENTIFICATION OF STAR FORMING 
GALAXIES 

Both star formation and AGN activity cause the host 
galaxy to (re-)emit at (mid-)IR and radio wavelengths. 
To study the IR/radio properties of these two distinct 
populations separately, information from different re- 
gions of the electromagnetic spectrum is thus required. 



ISmolcic et al.l (|2008l ) devised a method which, in a statis- 
tical sense, is capable of selecting star forming and AGN 
galaxies with a simple cut in re st frame optical col our. 
It relies on the tight correlation (Smolcic ct al. 200j|) be- 
tween the rest frame colours of emis sion line galaxies and 
their position in the BPT diagram (jBaldwin et al.lll981| ) 
and was developed and calibrated with radio sources at 
z < 1.3 using the principal component colour 19 hence- 
forth referred to as 'PI'. It can, however, be easily 
adapted to other rest frame colours because galaxy SEDs 
from th e near-UV to the NIR represent a one- parameter 
family (jObric et alJl2006t Islnofcic et al J 12006). Here we 
use the combination of the filters u and K to select AGN 
and SFGs. This choice is motivated by the desire to ap- 
ply the classification to both the radio- and IR-selected 
sample; the likely presence of dust-obscured star forming 
systems in the IR-selected sample requires the inclusion 
of a red band, to prevent, as best possible, dust-reddened 
star forming sources from being mistaken for red, early- 
type AGN host galaxies. 

In Fig. [71 we show the correlat ion of PI (computed ac- 
cording to ISmolcic et all 120081 ) and (it - K) for -950 
VLA-COSMOS sources, for which both PI and (u - K) 
were available, rest frame (u — K) colours were com- 
puted with ZE BRA (Zurich Extraqal actic Bayesian Red- 
shift Analyzer, Fel dmann et al.l [20061 ) which was used to 
find the best-fitting SED template to the COSMOS pho- 
tometry in the medium and broad band filters u* , B, 
V, g + , r + , i + , z + , J and K s , as well as in th e fou r 
IRAC channels given the known redshift (cf. § IA.5j) . 
Note that the magnitudes u and K used here are com- 
puted in Johnson-Kron-Cousins filters rather than the 
COSMOS filters. An ordin ary least squares (OLS) bi- 
sector fit (Iso be et al.l [1990) accounts for the fact that 
both colours are subject to uncertainty and returned a 
best-fit correlation given by 

PI = -0.94(±0.0006) + 0.45(±0.01) X (u - K ) , (3) 

which is indicated in g r ey in Fig. [7j The criterion PI < 
0.15 of ISmolcic et al.1 (120081 ) for the separation of SF 
(PI < 0.15) and AGN sources (PI > 0.15) thus cor- 
responds to (u — K) = 2.42. Note that due to our treat- 
ment of composite SF/AGN sources we adopt a slightly 
different colour threshold for the selection of SFGs (see 
following parag raph and Fig. [ 8 l) . 

From Fig. 24 of lSmolcic et a l. (2008) (reproduced in the 
upper left corner of Fig. [5]) it is obvious that the tails of 
the distribution of star forming and AGN systems in PI 
colour space extend into the selection regions for AGN 
and star forming sources, respectively. Moreover, about 
10% of the sample on which the classification scheme 
was developed are 'composite' systems and found on ei- 
ther side of the divide 20 . When a source is classified as 
SF or AGN based purely on its rest frame optical colour 
there thus is a non-negligible probability of assigning it 
to the false category. For some purposes, e.g. when 

19 PI and its homologue P2 are linear combinations of the 
narrow band (modified) Stromgrcn filter magnitudes (uz, vz, bz, 
vz; IQdell et al.l H2002I )) in the wavelength range 3500-5800 A; see 
Smolcic ct al. (2008) for the definitions and additional details. 

20 The cl assification of sources in the reference sample of 
ISmolcic et al.l l|2008T ) into AGN, SF, and composite galaxies is 
based on their positi on in the opt ical spectroscopic diagnostic 
(BPT) diagram ^Baldwin et al.|[T98ll) . 
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estimating which fraction of AGN systems have similar 
IR/radio properties as star formers, it is thus useful to 
adopt a probabilistic approach. Given the distributions 
^SF, -^Vagn and A compo . (see Fig. [8]a) a possible defi- 
nition for an effective probability Pr e ff. (SF) of correctly 
classifying a source as star forming at a given rest frame 
optical colour is 



Pr cff . (SF) « Pr (SF) + 



SF 



SF 



where 



Pr (SF) = 
Pr (compo.) : 



^Vagn 



N SF 



x Pr (compo.) , (4) 



N SF + N AGN + A compo . 
N 

1 v compo. 

iV S F + -/Vagn + ^V compo . 



In setting up (j4|) we have assigned composite systems to 
the SF and AGN population according to the relative 
abundance Asf/Aagn of SF and AGN sources at the 
particular colour. Pr e g . (SF) is a priori given as a proba- 
bility as a function of PI through the distributions Asf , 
-^Vagn and A compo . presented in ISmolcic et al.l {2008). 
However, it may be directly converted to the desired de- 
pendency on (u — K) by convolving the expression in Q 
with the distribution of PI at fixed (u — K) colour (see 
inset of Fig. [7]), which reflects the range of probabilities 
Pr cff . (SF | PI) that contribute to Pr eff . (SF | (u - K)). In 
the upper right panel of Fig. [8] we show the distribu- 
tions Pr c g. (SF/AGN|P1) obtained according to equa- 
tion Q and smoothed with a three-point running av- 
erage (black curve - star forming sources; dark grey 
curve - AGN systems). Its convolution with a stan- 
dard normal curve leads to the probability distribution 
Pr efI . (SF/AGN| (u - K)) shown in the lower panel of 
Fig. [H] which uses the same colour scheme as in panel 
(b). The uncertainty in the best-fit correlation between 
PI and (u — K) has been translated into an error in the 
probability function which is shown as a light grey area 
to either side of the black line giving Pr c ff. (SF | (u — Kj) 
in panel (c). Due to the small uncertainties in the 
OLS bisector line parameters of equation {3J the disper- 
sion o~api is the most important factor that determines 
the differences in the shape of Pr c ff (SF/AGN | PI) and 
Pr off .(SF/AGN|( U -X)). 

If one assigns composite objects to the SF and AGN pop- 
ulation according to equation Q the point of equal prob- 
ability of correctly classifying objects as SF or AGN, re- 
spectively, is reached at (u—K) = 2.36. This value is only 
slightly different from the direct translati on (see previ- 
ous pa ragraph) of the original definition in ISmolcic et al.l 
( 2008). In the remainder of the paper we will use the 
(u — K) = 2.36 threshold to separate SFGs from sources 
with emission that is dominated by AGN activity 21 . 
Apart from the tails in the colour distribution of AGN 
and SF systems which cross the colour threshold, three 
additional effects could reduce the accuracy of the clas- 
sification scheme. 

First of all, a general evolution of the SF and AGN pop- 
ulation to bluer colours at high redshift would lead to in- 
creasing contamination by AGN of the high-z population 

21 When writing about and plotting probabilities we will hence- 
forth use Pr (SF) as a shorthand for Pr eff . (SF | (u - K)). 



of SFG s if the colour cut is not adapted. ISmolcic et al.l 
(2008) have shown that an unchanging threshold is ad- 
equate until at least z ~ 1.3. In Fig. [9] we plot the 
distribution of (it — K) colours of our sources and follow 
the evolution of the median colour for both the IR- (dark 
grey histogram) and the radio-selected sample (light grey 
histogram; sources common to both samples lie within 
the hatched area) . We find no evidence for a strong evo- 
lution of average colours in either of the two samples out 
to z ~ 2, and out to z ~ 3 only by a small amount. Hence 
we apply the selection criterion uniformly to all sources, 
regardless of their redshift, except for the objects at the 
highest redshift where the medians have begun to change 
appreciably (see lower right panel of Fig. [9]). 
Secondly, nonperiodic flux variations of active galaxies 
will affect the choice of the best-fitting SED if photo- 
metric measurements are not simultaneously carried out 
over the whole spectrum. Since the rest frame optical 
colours are determined using SED templates this can 
cause misclassifications of AGN or SFGs with colours 
close to the threshold (u — K) = 2.36. A variability anal- 
ysis (M. Salvato, private communication) of our 1.4 GHz 
sources revealed that maximally 20% of these display 
stron g variability (defined as VAR > 0.25; cf. equation 
(1) in Sal vato et al.ll2009ft . The true fraction of affected 
sources is likely to be smaller because inaccuracies in the 
photometry can artificially raise the variability parame- 
ter. 

Finally, we cannot exclude that some unobscured Type 1 
AGN with a blue (u—K) colour will be assigned to the SF 
category in our classification scheme. It is also possible 
that a number of dust-reddened starburst galaxies end 
up being classified as AGN, even though we used a red 
filter to define our rest frame colour on which we base 
the separ ation into SFGs and AGN. In the calibration 
sample of Sm olcic et al.l (|2008ft this kind of contamina- 
tion amounted to less than ~10% (see their Appendix 
B2). 

4. IR SED TEMPLATE FITTING 



Data from lensed high-z galaxies (|Siana "et~all[200l 
IGonzalez~e t al. 2009]) and from recent deep FIR surveys 
have shown that the SEDs of local star forming galax- 
ies reproduce the SE Ds of high redshi f t galaxies well 
out to z < 1.5 (e.g. lElbaz et al.l l2002t iMagnelli etahl 
120091: IMurphv et al.ll2009b[ ). However, it has also been 
reported that the SEDs of some IR-selected galaxies 
at high redshift can differ from loca l templates both 
at M IR (jRigbv et al.l 12008ft and FIR (|Svmeonidis et al.l 
12008ft wavelengths, conceivably due to intrinsic scatter 
in the physical properties of these sources which deviate 
from the median trend that the empirical galaxy tem- 
plates represent. 

Followi ng the procedure described in IMurphv et al.1 
(|2009bft . we derive infrared luminositie s (Z>tir) by fit- 
ting th e 24 and 70 um data points to the lCharv fc Elba3 
(|2001ft SED templates and integrating between 8 and 
1000 /im. This wavelength range is in principle also 
sampled by S-COSMOS observations a t 8 and 160 ^m 
(|Sanders et al.ll2007t IPraver et al.l [2009) but we restrict 
ourselves to the two aforementioned bands because (?) 
at z > 0.6 a 8 am measurement would include stellar 
light (which starts to dominate at the SED at rest frame 
5 /mi), while we are fitting pure dust templates; and (it) 
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the shallower coverage and broad PSF of the 160 /xm ob- 
servations complicate the identifi cation of unamb i guous 
counterparts. Our choice of the iCharv fc Elbazl (2001) 
templates is motivated by the fact that they have been 
found to exhibit 24 /70 nm flux density r atios that are 
more representative (jMagnelli et al.ll2009t ) of those mea- 
sured for galaxies at z ~ 1 comp ared to the lDale fc Heloul 
(|2002[ ) or lLagache et all (|2003f) templates. 
For the cases where a source is detected firmly at 24 
and 70 ^tm, the best-fit SEDs are determined by a \ 2 
minimization procedure whereby the SED templates are 
allowed to scale such that they are being fitted for lu- 
minosity and temperature separately. Consequently, the 
amplitude and shape of the SEDs scale independently to 
best match the observations. The input photometry is 
weighted by the S/N ratio of the detection if it is a well- 
defined measurement, and the normalization constant is 
determined by a weighted sum of observed-to-template 
flux density ratios for all input data used in the fitting. 
In the cases where only an upper limit is available at 
70 /im, the latter is not incorporated into the x 2 mini- 
mization but used to reject fits which have flux densities 
above the associated measured limit. 
Errors on the best-fitting value of Ltir are determined 
by a standard Monte Carlo approach using the photomet- 
ric uncertainties of the input flux den sities which reflect 
both calibration errors (~2 % at 24 (lEngelbracht et al.l 
[2001 and -5% at 70 /mi ([Gordon et alj 120071 )) and 
the uncertainties in the PSF-fitting (generally of order 
Spsf/(S/N), where S PS f is the flux returned by the PSF 
fit). 

If a source is only detected at 24 /im, we also fit the 
photo metry using the SED templates of iDale fc Heloul 
( 2002) and define the best estimate of the IR luminosity 
as the average Ltir from the two separate fits. 

5. SELECTION EFFECTS AND STATISTICAL 
TREATMENT OF FLUX LIMITS 

5.1. Shifts Between the Average IR/Radio Ratios of 
Flux Limited Samples 

The selection effects that are the topic of this section 
arise in flux limited samples when flux information from 
one of the selection bands is directly used in the com- 
putation of the quantity being studied. In the present 
case the critical quantity is the logarithmic IR/radio flux 
ratio </, but analogous effects need to be considered in 
the context of studies of the distribut ion of spectral in- 
dices at differe nt radio frequencies fe.g. lKellermannlll964l : 
lCondonlll984D of X -rav to optical continuum slopes of 
AGN (lFrancislll993l) or of th e M. - a and M. - L rela- 
tionships ([Lauer et al.ll2007j) . 

In Fig. (TU] we illustrate the origin of the selection ef- 
fect: consider the left hand panel in which the IR-to-radio 
SEDs of three sources with different observed bolomet- 
ric flux are distributed along the vertical axis. Each of 
these three SEDs splits into three branches at the peak 
of the SED, thereby schematically reflecting the range 
of observed IR/radio ratios (from top to bottom: 3 a 
radio-excess outlier - dashed line; average source - solid 
line; and 3 a IR-excess outlier - dotted line). If we im- 
pose the indicated selection threshold at 1.4 GHz (red 
line) the resulting sample will contain (?) all sources of 
the brightest flux class, regardless of their IR/radio ra- 



tio; (ii) the source with an average IR/radio ratio and 
the radio-excess source from objects of the intermedi- 
ate flux class and; (Hi) in the faintest flux bin only the 
radio-excess sources. Since the fainter sources are more 
abundant (as parametrized by the slope of the differen- 
tial source counts dN /dS oc S~@, with (3 > 0) this results 
in a surplus of radio-excess sources and consequently a 
low average IR/radio ratio in a radio-selected sample. 
The right hand side of Fig. [TU] shows that an IR-selected 
sample is biased in the opposite direction, i.e. towards 
high IR/radio ratios. 

The analytical expression for the difference between the 
average IR/radio ratio of IR- and radio- selected sam- 
ples is (iKellermannl H96l iCondonl H98l iFrancisI lilM 
ILauer et al.ll2007t) : 

Ag bias = ln(10) (fi - 1) a\ . (5) 

It thus depends on /?, the power law index of the source 
counts, and on a qi which is the dispersion of the IR/radio 
relation. Note that this offset will occur regardless of the 
relative depth of the two involved bands. An estimate 
of the 'intrinsic' (i.e. unbiased) IR/radio ratio can be 
obtained by constructing the sample using an unrelated 
selection criteri on like optical lum inosity, mass or mor- 
phological type (jLauer et al.ll2007l) . 
Since the recent work on the evolution of the IR-radio re- 
lation at intermediate and high redshift was often based 
on flux limited surveys, we would expect most of the 
findings to be affected by this selection bias to a certain 
extent. In Table [5] we have collected the selection cri- 
teria and average values of q (final column) that were 
published in the literature during the last decade. 
We see that broadly speaking the various IR/radio diag- 
nostics have values </24 ~ 1, 970 ~ 2.1, <7fir. ~ 2.3 and 
<7tir ~ 2.6. These different values are not the result of 
selection effects but reflect if the IR filter covers a wave- 
length range that is close to the IR SED peak or a part 
of the SED with lower energy content. In the following 
paragraph we will discuss the plausible influence of selec- 
tion effects on the various measurements of (724 and (/fir, 
in particular. 

Due to the high sensitivity of the MIPS 24 /im band many 
of the papers listed in Table [5] have studi ed the IR/radio 
ratio </24- The radio-selecte d samples of lAppleton et al.1 
( 2004) and llbar et all (gOOS) find that g 2 4 € [0.94, 1], de- 
pending on the choice of the IR template use d for the K- 
correction. The local IR-selected sample of IRieke et al.1 
(2009) on the other hand has a mean ^24 of ~1.25 
and shows some signs of variations with IR luminosity. 
The offset between the means of the radio-selected sam- 
ples and the IR-selected data set is ^0.3 dex, in good 
agreement with the predicted A^bias = 0.31 of equa- 
tion ([5]) if we set a a sa 0.3 in accordance with obser- 
vatio ns (e.g. lYun et al.l l2001h [Belli 120031: lAppleton et~aH 
12004 llbar et al.ll2008fl and under the simplified assump- 
tion of Euclidean source counts ((3 — 2.5). In the 
case of the FIR/radio flux ratio (/fir we can c o mpare 
the two radio-sel e cted ( sub)samples of IGarrett] (2002) 
and ISaiina et al.l (|2008f ) that have </fir ~ 2 with a 
jointly radio- and s ub-mm selected mean of 2.07 from 
IKovacs et al.l (120061) and mean values qpm. 6 [2.2, 2.4] 
for IR-selected ([Younger et Id] [2001 ISaiina et al.ll2 008) 
or essentially volume limited samples in lBelll (|2003|) and 
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lYun et all (|200I| ) . As with g 2 4 there is thus evidence of a 
~0.3 dex shift in (?tir between radio-selected and other 
samples. As far as we know no measurement of qro in 
an IR-selected sample exists but the compilation in Ta- 
ble shows that reassuringly all determinations of gyp 
based on radio-se l ected samples ( Appleton et al.l 120041 : 
iFraver et alll200l iSevmour et all 12009) are quite simi- 
lar. 

The radio sta c king exper iments of iBovle et all (120071 ) , 
iBeswick et all (|2008f ) and iGarn fc Alexander! (|2009D do 
not fit the picture which is probably due to the differ- 
ent nature of the analysis. Nevertheless, it cannot be 
excluded that part of the variations in the other studies 
are due to field-to-field variance or different assumptions 
about IR SEDs and the radio spectral slope. To this 
end we will test in §|6] whether or not the offset between 
our IR- and radio-selected samples - that have been con- 
sistently constructed from the same parent data sets - 
conforms to our expectation. If so, it would be strong 
support for selection effects alone being able to recon- 
cile the seemingly discrepant measurements of average 
IR/radio properties in the literature. 

5.2. Derivation of Distribution Functions with 
Survival Analysis 

Discarding the information from undetected counter- 
parts introduces a second source of bias in addition to 
the selection effects mentioned in § 15.11 It arises from 
the unrepresentative sampling of the true distribution 
function of IR/radio ratios by sources which are directly 
detected in all involved bands. We would like to empha- 
size that the shift Aqbias in equation §5§ is the difference 
between the mean q of IR- and radio-selected samples 
with correctly sampled distribution functions. Agbias 
cannot be compensated by accounting for upper or lower 
limits on q due to undetected IR or radio counterparts 
in the two different samples; as discussed in the previous 
section, the mean IR/radio ratio measured in an IR- 
and radio-selected sample only brackets the value one 
would measure with an unbiased data set which we can 
best approximate by a sample jointly selected at IR and 
radio wavelengths (see § 12.3. 3|) . 

The ratio of two flux constraints that could be either a 
well-defined measurement or an upper limit will render 
an upper or lower bound, a well-defined value or be inde- 
terminate (if both numerator and denominator are lim- 
its). Since we use the pooled information from a radio- 
and IR-selected sample in this study the latter case never 
occurs. We do expect, however, to encounter upper lim- 
its on IR/radio ratios from radio-selected sources that 
are not detected in the IR or lower limits if the radio 
counterpart of an IR - select ed source was too faint to be 
detected (cf. SSlA~T1fclA3|). 

Let qi (i = 1, ... , n) be the actual value of the flux ratio 
for each of the n sources in a suitably defined sample 
(e.g. the population in a certain slice of redshift). As a 
consequence of the noise characteristics in the radio and 
IR images qi can only be measured if it lies in the interval 
[<h i o£]i where q~ and qf are upper and lower limits on 
the flux ratio, respectively. These limits may be different 
for each source. Our knowledge about the distribution of 
IR/radio ratios after carrying out all our measurements 
can thus be summarized with two vectors of variables, Q 



and 8: 



Qi= max(min((ji, qf), q i ) 



(6) 



-I \iqi<q[ 

if ft e [qT,Qi] 

1 if <?i > qf 



In survival or life time analysis the action of imposing 
measurement limits is referred to as 'censoring'. A vari- 
able Qi is said to be left censored if Qi < q~ and right 
censored if Qi > qf . If both kinds of censoring occur in a 
data set it is called doubly censored, otherwise one talks 
of single censoring. During the remainder of the paper 
we will use the terms 'limit' and 'censored measurement' 
interchangeably. 

In Appendix [C] we sketch the steps that are involved in 
going from the information (Q, 8) to the distribution 
function of the IR/radio ratios. Inferring the true distri- 
bution of the qi of a sample is essential for the calculation 
of its average IR/radio properties. In §|6] we will con- 
struct distribution functions for data sets that are both 
singly and doubly censored. Recipes for dealing with 
the former case are plentiful in texts on survival analy- 
sis (see e.g. iFeigelson fc Nelson! |l985) for applications to 
astronomy) such that we only include some brief remarks 
in § IC.2I Since the more general case of double censor- 
ing is not as widely used in astronomical applications, 
the most important formulae and u seful computational 
guidelines are provided in Appendix lC.fi 
The methods described in Appendix [C] have been im- 
plemented using Perl/PDL 22 scripts written by M.T.S. 
. Their correct functionality was verified with examples 
in the literature. In particular, we checked that our im- 
plementation of the algorithm for the calculation of t he 
doubly censored distribution function (|Schmitt]|1985j) - 
when applied to the special case of singly censored data - 
gave the same results as the scri pts based on the Kaplan - 
Meier product lim it es timator (|Kaplan fc Meierl (|1958|) : 
see also Appendix IC.2|) . 

6. RESULTS 

The main focus of this section is the search for changes 
with redshift of the average IR/radio ratio in the SF pop- 
ulation. We track evolutionary trends in the range z < 
1.4 for both monochromatic and TIR/radio flux ratios in 
§§ m and !6.2l and separately consider a sample of high 
redshift (z > 2.5) sources in § 16.51 § 16.31 is dedicated to 
the IR/radio properties of AGN hosts and in § 16.41 we 
study variations of IR/radio ratios with luminosity. 
Previous studies have carried out similar analyses us- 
ing a variety of IR/radio diagnostics. These include 
MlPff-based monochr o matic flux r atios (? 24 and gyp (e.g. 
lAppleton et al.l (|2004D ; IIbar et all ([20081 ) ; ISevmour et all 
(2009) ; see equation ^ for the definition of q 2 4 po] ) which 
we discuss in § 16.11 Other studies have used the FIR 
(42 . 5-1 22.5 pn) to radio flux ratio of t r. (see equation (pi ; 
e.g. lGarrettll2002l:IKovacs et al.ll200alSaiina et al]l2008f) . 

22 The Perl Data Language (PDL) has been developed by K. 
Glazcbrook, J. Brinchmann, J. Cerney, C. DeForest, D. Hunt, T. 
Jenness, T. Luka, R. Schwebel, and C. Soeller and can be obtained 
from http:/ /pdl. perl. org 
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or the ratio of total IR luminosit y (Ltir) to radio lumi- 
nosity (e.g. |Mur^^ir23llQM9): 

gTIR = log ( 3.75^0" w ) ' log (wip) - (7) 

Total infrared luminosities Ltir (in units of [W]) 
are calculated by integrating the SED between 8 and 
1000 /im. The rest frame 1.4 GHz luminosity (expressed 
in [WHz- 1 ]) is 

Liagb* = f^ L [f a S^lAGRz) , (8) 
(1 + z) L 

where S V (1A GHz) is the integrated radio flux density of 
the source and Dl(z) the luminosity distance. The K- 

correction Kiaghz(z) = (1 + z)^ 1- ") depends on the 
spectral index a of the synchrotron power law S v oc v~ a . 
For the rest of the analysis we will assume that a = 0.8 
(|Condonl lT992). We will return to the TIR/radio flux 
ratios in § 16.21 

6.1. Monochromatic IR/Radio Properties of Star 
Forming and AGN Galaxies 

6.1.1. Observed Flux Ratios 

The observed 24/xm/1.4GHz flux ratio (724 obs is plot- 
ted against redshift in Fig. [ITI for SFGs (top) and AGN 
(bottom). Sources are assigned to the two categories de- 
pending on whether Pr (SF) is larger or smaller than 50% 
(cf. ®. 

While there clearly are many radio-loud sources in our 
AGN sample, Fig. [TTJb shows that a majority of the ob- 
jects assigned to the AGN category displays very similar 
IR/radio ratios as the SFGs. We will discuss this ob- 
servation in more detail in § 16.31 At the same time, the 
sample of SFGs also includes a number of radio-excess 
sources. They usually have photometric redshift esti- 
mates and mostly lie at 1 < z < 3. This roughly corre- 
sponds to the redshift range in which photometrically de- 
termined redshifts are subject to the largest uncertainty 
because the 4000 A break is only sampled by broad and 
widely spaced photometric bands. As a consequence, ab- 
sorption features and emission lines from AGN and SF 
systems often interfere with each other in the same fil- 
ter. Even though we did attempt to remove all unreliable 
redshifts - as described in 5 IA.5I - it thus seems proba- 
ble that at least some of these cases are due to wrong 
distance estimates and hence to the selection of an in- 
appropriate optical SED. Since this results in a faulty 
(u — K) colour, the source in question could then have 
been assigned to the SF rather than the AGN category. 
Another po ssibility is that the peak of AGN activ ity at 
z ~ 2 (e.g. IWolf et al.l 120031; IRichards et all [2001 also 
influences the SF sample due to the statistical nature of 
the identification of SFGs and because especially AGN 
in composite systems could have been classified as star 
forming. Lastly, we tried to assess if unobscured Type 
1 AGN represent a significant fraction of the nominally 
star forming radio excess sources in the perti nent r cd- 
shift range. Based on the confidence class (see 5 IA.5.1I) of 
those objects for which follow-up spectroscopy was avail- 
able, we estimate that only ~5% are quasars classified as 
SF due to their blue colour. 



The mean value of 924. obs decreases as a function of red- 
shift. We will show later on (see Fig. [20]) that this de- 
crease agrees with the variations local LIRGs (detectable 
only out to z ~ 1; see vertical dotted lines in panel (a)) 
and ULIRGs would display if redshifted. 
In Fig. [12] we plot the observed 70/xm/1.4GHz flux ra- 
tio 970. obs of our sources as a function of redshift. All 
symbols and colours are exactly as in the previous fig- 
ure. Note that censored measurements due to flux limits 
at 70/mi are more frequent than was the case for 924. obs 
because the 70 ^m observations are much shallower. The 
observed flux ratio before JT-correction shows the same 
decline at higher redshifts as was seen for the observed 
24/jm/1.4GHz flux ratios. 

6.1.2. Evolution of 524 and 570 with Redshift 

As described in §[4] all sources classified as star 
forming were fit with IR SEDs in order to derive IR 
luminosities. As a by-product of the template fitting we 
can immediately obtain rest frame (i.e. -fT-corrected) 24 
and 70 fim flux densities by convolving the best-fitting 
SED with the filter response functions of MIPS. In 
the following we define the if-corrected 24 and 70 /im 
fluxes as th e average of the value s obta ined from the 
librari es of iCharv fc Elbazl (f200lh and iDale fc Heloul 
(2002) and use them to construct monochromatic rest 
frame flux density ratios (724.0 an d <?7o. o- The associated 
1.4 GHz flux densities have been ^-corrected according 
to equation ((8|). 

To quantify the evolution of c/24, in the joint IR- and 
radio-selected sample we: 

1. Bin the data such that each redshift slice contains 
an identical number of objects (~250). The num- 
ber of bins is kept limited to guarantee that the dis- 
tribution function of (724, is sampled sufficiently. 

2. Run the iterative algorithm outlined in § IC.1I to 
find the cumulative distribution function of (724, 
at each redshift. The median immediately fol- 
lows from this computation as does the scatter in 
the population which we obtain by fitting a Gaus- 
sian distribution with known mean (equal to the 
previously determined median) to the distribution 
function. The choice of the Gaussian is moti- 
vated by the shap e of the local IR-radio relation 
IjHelou et al.lfl985l : lYun et al.ll200TI : IBeflll2003l) . 

3. Determine the evolution of the average IR/radio 
ratio by fitting a linear trend line to the medians. 
Only measurements at z < 1.4 are considered for 
this since the scatter at higher redshifts is found to 
increase abruptly, thus making the determination 
of the median uncertain. 

Steps 1 and 2 are are also carried out individually for 
the sample of IR- and radio-selected galaxies. The re- 
sults are shown in Fig. 1131 Since the cumulative distri- 
bution function is normalized it lies in the range between 
zero and unity and can thus be regarded as the proba- 
bility of obtaining a measurement of (724, which is less 
than - in the case of the radio-selected sample (light grey 
curve) - or in excess of - for the IR-selected sample (dark 
grey curve) - the ordinate. The distribution function of 
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the doubly censored union of the IR- and radio-selected 
samples is also parametrized such that it runs from to 
1 with increasing 924. o- It is plotted in black together 
with a dashed red line which shows the corresponding 
best-fitting Gaussian distribution. The intersection of 
the black curve with the 50% probability line (dotted 
horizontal line) defines the median value of 924. • 
Fig. [13] demonstrates that the median of the radio- 
selected population lies systematically below that of the 
IR-selected objects. The shift is approximately 0.35 dex 
at low redshift and grows to about 0.7 dex beyond z ~ 
1. The increase is probably caused by the intrinsically 
high er scatter in the I R-radio relation at h igh luminosi- 
ties (|Yun et all [20011: iBressan et all 120021 ). possibly in 
combination with the reduced reliability of photomet- 
ric redshifts and/or some falsely classified AGN that be- 
gin to affect the sample starting at z ~ 1.3. A shift of 
~0.35 dex as observed at z < 1 where the accuracy of 
our measurements is highest agrees fairly well with the 
prediction of equation ((5J and is hence a likely explana- 
tion for differences between previously reported average 
IR/radio properties of both local and high-z galaxies (e.g . 
lAppleton et~aLll2004t llbar et al.ll2008b iRieke et al.ll2009f ). 
In the lower panel of Fig. [13] we plot the medians (924, 0} 
of the jointly IR- and radio-selected SFGs (black dots) at 
different redshift on top of the if-corrected values (724, 
(colours and symbols are identical to those in Fig. [TTIa). 
The error bars mark the 95% confidence interval associ- 
ated with the median. Table [3] lists the median and scat- 
ter of (724, which were determined with survival analysis 
in each of the redshift bins of Fig. [T31 In addition to the 
measurements carried out on the jointly IR- and radio- 
selected sample the table also provides the according val- 
ues for the IR- and radio-selected samples individually. 
Using the 2 a errors on the medians as weights we fit 
them with a model of linear redshift evolution. The best- 
fitting trend of ((724, 0) y S- z is shown in black in the lower 
window of Fig. [T3](see Table@]for the parametrization of 
the line). Because the fit was carried out with respect to 
linear redshift space while the plot has a logarithmically 
scaled redshift axis it is curved. Within the errors the 
slope d (i24,a)/dz = —0.015 ± 0.136 is consistent with no 
evolution of the i-T-corrected 24 /zm/ 1.4 GHz ratio at z < 
1.4 (the maximal distance out to which the precision of 
the photometric redshifts is high). The y-axis intercept 
of the tren d line at z = is in agreement with the recent 
analysis of Ri eke et al.l (|2009l ) who find that (724,0 = 1-22 
with a scatter of 0.24. In our sample we find ((724, o)z=o — 
1.28±0.10 (where the error states the formal la uncer- 
tainty from the linear fit). For a comparison between 
the average IR/radio properties of radio-selected sam - 
ples we can refer to th e studies of lAppleton et all (]2004f ) 
and llbar et al.1 (|2008f l: they report an average 524, — 
0.94- 1, depending on the IR-SED adopted to if-correct 
to the rest frame. These values agree well with the range 
of medians ((724,0) £ [0-8, 1] measured for radio-selected 
COSMOS data at intermediate redshift (cf. left-most 
column of Table [3]) . 

Our convention for choosing SFGs states that Pr (SF) 
must be at least 50%. In Fig. [14] we assess how chang- 
ing this threshold affects the redshift evolution that is 
inferred from the data. The variation of the parame- 
ters of the best-fitting evolutionary trend line is shown 
in the upper- and lower-most window of Fig. [14] (y-axis 



intercept and slope, respectively) . A black symbol in the 
middle of the displayed data range marks the results that 
were shown in Fig. 1131 They are fully consistent with the 
evolution found if a more conservative threshold - e.g. at 
Pr (SF) = 66% - for the selection had been chosen. It is 
interesting that the inclusion of a significant fraction of 
sources with a probability of up to 80% of being AGN 
does not alter the results either. This is a strong indica- 
tion that our sample of optically selected AGN contains 
many objects with IR/radio properties that closely re- 
semble those of SF systems. Sim i lar ob serva tions were 
made by, e.g.. iSopp fc Alexanderl (|1991l ) and IRov et al.l 
(1998), who studied local samples of radio-quiet quasars 
and/or Seyfert 1 sources lacking a compact nucleus. The 
middle row of Fig. [i"4l shows that while the average val- 
ues of (724, are similar for many SFGs and AGN the 
latter are sub ject to a large r scatte r as was pre viously 
found by, e.g.. ICondon et al.1 (fl98l : lObric et al.l pOOfo : 
IMauch fc Sadlerl (120071) . ~ 

As an additional test of the robustness of our findings we 
checked if the evolutionary trend in SF samples selected 
through (u— K) or PI differs. For the radio-selected sam- 
ple where both colours were available we found equivalent 
results regardless of the chosen approach. 
Figures [T5] and [TO] (the results of which are summarized 
in Table [5]) repeat the analysis of Figs. [T3] and [TJ] for 
the if-corrected 70 /im/1.4 GHz flux ratio. Note that in 
comparison with Fig. [12] the number of censored mea- 
surements is much smaller when we consider rest frame 
IR/radio ratios 970,0 rather than observed flux ratios 
970. obs- The reason is that we do not require a 70 /im 
detection for the IR template fitting but also fit objects 
which have a limit at 70 /im and a direct detection at 
24 ztm. 

The plot of (770, vs - redshift (lower panel of Fig. [15]) 
as well as Fig. [TBI which illustrates the stability of 
the findings with respect to changes in the selection cri- 
terion for SFGs, show that (770, behaves in a similar 
way as was found for 924,0- The extrapolated average 
70/xm/1.4GHz flux ratio at z = 0, ((770,0)2=0 equals 
2.31±0.09 and if AGN are included the scatter in the rela- 
tion increases in analogy to what was found for 924, o- As 
for q 2 4, the evolutionary slope d (i™,o)/dz = — 0.123±0.135 
(slope and normalization of the evolutionary trend line 
are logged in Table 0]) is consistent with zero. 
The average IR/radio prop er ties o f the r adio-selected 
sampl es of lAppleton et al.l (|2004l ) and iFrayer et al.l 
([2006]) 23 are 2.16±0.17 and 2.10±0.16, respectively. Al- 
though the agreement with our findings is not quite as 
good as in the case of 924, they are are nevertheless con- 
sistent (within both the formal error and the scatter) 
with the range of medians (970,0) £ [1-7, 2.1] at interme- 
diate redshift in our radio-selected sample (see Table [5]) . 
To our knowledge there so far has been no comparable 
study which uses an IR-selected sample to compute an 
average 970. 

6.2. Evolution of TIR/Radio Flux Ratios with Redshift 

23 Although based o n a catalog of xFLS 70 /an sources the sample 
of Fraycr et al. (2006) becomes essentially a radio-selected sample 
at the stage when sources with out a counterpart in the 1.4 GHz 
radio catalog of the FLS field (ICondon et al.1 2003) are removed 
from the sample. 
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In the local universe the correlation of IR and radio 
flux is tightest if integrated (F)IR luminosities rather 
than monochromatic flux ratios are considered. To com- 
plement the analysis of § 16. II we thus show in this section 
the correlation of TIR (8-1000 fixn) and 1.4 GHz luminos- 
ity as parametrized by the TIR/radio ratio 9tir for our 
VLA- and S-COSMOS data. 

The computation of the distribution functions for the 
parameter 9ti r is c arried out following the same steps 
described in § 16.1.21 The results are shown for a num- 
ber of redshift bins in Fig. [17] where we also compare 
the median derived for the jointly IR- and radio-sel ected 
SFGs with the local value of (9tir)z=o = 2.64±0.02 pell 
(2003); vertical dashed line). Our average values (qtir) 
in the range z < 1.4 (see Table [6]) lie to either side and 
always rema in well with in the dispersion of the local mea- 
surement of Eel (|2003| ). 

The evolution of 9tir is shown in the lower panel of 
Fig. H7] using the same presentation of the data as for 
the monochromatic IR/radio flux ratios. Since the lat- 
ter were derived based on the IR templates which are 
used here to calculate the integrated IR luminosity, we 
expect by construction that the evolutionary trend is in 
good qualitative agreement with the findings of § 16.11 For 
the same reason we cannot expect to observe a reduced 
scatter in the values of 9tir with respect to those of the 
monochromatic flux ratios as the spread in the properties 
of the best-fitting IR SEDs must manifest itself in (724, 
and <77o, as well. 

The line parameters for the evolution of (<?tir) are given 
together with those of (924,0) and (qro,o) in Table 01 
in contrast to 524,0 and #70,0 the best-fitting evolution- 
ary trend for (<7tir) suggests a decrease of the average 
TIR/radio ratio by 0.35 dex out to z ~ 1.4. However, 
this slope is detected at the 2 a significance level and pre- 
dicts a median (<?tir) at z ~ 1.4 that still lies within the 
dispersion measured in our lowest redshift bin. It thus 
seems unlikely that the evolutionary signal is real, espe- 
cially in view of the results of § !6.5l where we measure an 
average (9tir)z>2.5 that is in excellent agreement with 
the local value for a subset of highly redshifted galax- 
ies in the COSMOS field. An examination of the evo- 
lutionary slopes for (924,0), (970,0) and (qtir) in Table 
0] shows that they become more negative along this se- 
quence. This could be related to a number of radio-excess 
sources with 1 < z < 3 which are part of our optically se- 
lected SF sample (visible as a diffuse cloud of upper limits 
and detections below the main locus of symbols in all our 
plots of q vs. z; see also our comment in § I6.1.1[) and that 
tend to lower the average IR/radio in this redshift range. 
If these objects were falsely classified composite sources 
or AGN the increased emission at 24 /im from their hot 
dust might be able to compensate the radio-excess, thus 
leading to zero evolution in 924, as observed. 970. and 
9tir on the other hand sample mainly IR light from star 
formation and hence are lowered in the presence of ex- 
cess radio emission. This scenario can also explain why 
the evolutionary slope of (924.0) is insensitive to the se- 
lection criterion for SFGs (cf. Fig. IT3)) while it varies in 
the same sense as described above in the case of (970, 0) 
and (9tir). 



6.3. AGN with Similar IR/Radio Properties as Star 
Forming Galaxies 

The analysis of the previous sections revealed (cf. Figs. 
HI ttH and US also Figs. QT] and [TJ that the IR/radio 
properties of SFGs are shared by many of the AGN- 
bearing systems in our sample. In this section we will 
study this in more detail. We first test f§ !6.3.1|) if it re- 
mains valid for a subsample of sources which are detected 
in X-rays and have been found to host an AGN using a 
different approach than the classification scheme intro- 
duced in §[21 In § 16.3.21 we then compute (in different 
redshift bins at z < 1-1.4) the relative frequency of AGN 
and SF sources as a function of the IR/radio ratio. 

6.3.1. IR/Radio Properties of X-ray Detections 

At the sensitivity of the XMM-Newton observations of 
the COSMOS field a large fraction of the detected sources 
is expected to be po wered by AGN. This is confirmed by 
ISalvato et~aT1 (|2009f ) who have shown that ~70% of the 
XMM-Newton sources have UV to NIR SEDs which con- 
tain an AGN component. In Fig. [19] we compare the 
observed 24 fim and 70 fim to radio flux density ratios 
924. obs and 970. obs of X-ray detected AGN hosts at dif- 
ferent redshifts with the predicted IR/radio properties of 
model SFGs (coloured tracks 24 ; see footnote and the text 
of § 16.3.21 for a dditional detai l s). No te that according to 
the analysis of ISalvato et al.l (|2009f) the AGN contribu- 
tion to the UV-NIR SED exceeds 50% for most of these 
sources. From Fig. [19] it is obvious that a majority of 
the XMM-Newton sources have IR/radio ratios that are 
perfectly consistent with those expected for starbursts. 
They are genuine examples of active galaxies in which 
the AGN, although significantly contributing to the SED 
at optical and X-ray wavelengths, does not cause signif- 
icant excess radio emission. Fig. [19] therefore is strong 
evidence that the findings of §§ 16. II and 16. 21 cannot be as- 
cribed to an inadequacy of the method we adopted to 
distinguish between AGN and SFGs. 

6.3.2. The Relative Abundance of AGN and SFGs on the 
Star Forming Locus 

We first define the 'main' locus of SFGs in a plot of 
observed IR/radio ratio vs. redshift. Working with ob- 
served flux densities is necessary because we want to 
avoid imposing a template fit with the IR SED of a 
SFG on an AGN-bearing source even if it quite prob- 
ably shares similar IR/radio properties. At each redshift 
the star forming locus is centred on the average value - 
(924 [70], template) - of 924, obs as predicted by the model 
SEDs of sources in the observable range of IR luminosi- 
ties. We then consider a region between +2 a q and -2 a q 

24 The tracks are constructed by taking the ratio of the K- 
corrections between (i) the flux density at the rest frame (A) 
and redshifted effective wavelength (A/(l + z)) of the MIPS fil- 
ter and, analogously, (ii) that applied to the 1.4 GHz band. The 
if-correction K(z) is defined as the ratio of the rest frame lumi- 
nosity Lj,(A) and the luminosity at wavelength A/(l + z) which is 
sampled by the observer's measurement of the flux ,S„(A): 

L„(A) = L„(Vi+») x K(z) = 4wD L (z) 2 S„(A) x K(z) . (9) 

Here D^(z) is the luminosity distance. The ^-corrections used in 
the conversion of observed IR flux measurements at 24 and 70 (im 
to rest frame quantities depends on the shape of the IR SED of the 
galaxies. For the radio flux it has the form given after equation 
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around (924 [70], template) in which we chart the relative 
frequency iagn (<Z, z) of AGN. The analysis is restricted 
to this ±2 a band because beyond it the sparse sampling 
of the distribution function of the IR/radio ratios leads 
to unwanted fluctuations of fAGN (<JS z) . The value oio~ q 
is a representative average of the scatter in our data for 
the SF population at z < 1, i.e. 0.35 dex. 
In the upper panel of Fig. [501 (Fig. HH shows the same 
information for q-jq) we show the expected variations in 
924. obs of high-z galaxies assuming that their SEDs at 
IR and radio wavelengths are similar to those of local 
SFGs. SEDs from three different template libraries - as 
well as that of the starburst M82 - are shown for differ- 
ent £tir. The tracks are normalized at z — using the 
best-fit evolutionary trend line displayed in Fig. 1131 In 
the background we re-plot (cf. Fig. [TTJa) the observed 
24 /im/1.4 GHz ratios of our sample of SF galaxies in or- 
der to show how they nicely follow the tracks of the local 
SEDs. The solid black lines delineate the ±2cr band cen- 
tred on ((724, template)- The jumps at z ~ 0.5 and 1 occur 
because the averaging of the IR templates is performed 
with a discrete and restricted set of IR templates. Since 
we merely use these boundaries to define the parameter 
space for the subsequent analysis the discontinuities are 
inconsequential. 

The expression for the relative AGN abundance which 
accounts for censored measurements and the use of dis- 
crete probability bins is (see derivation in Appendix |D|) 



fAGN (q,z) 



EI' 



A, 



i.etf. 



x (l-Pr(SF)),; 



(10) 



Here the summation with respect to i extends over a 
finite number n of probability bins. In the case of g24 
we grouped sources into bins of width APr (SF) = 0.1 in 
order to have a sufficient number of measurements, and 
thus to ensure a well-behaved estimate of the distribution 
func tion f(q) (computed according to equations (|Clj) and 
(|C2p ) in each probability bin. 

In the lower panel of Fig. [50] we present the function 

^AGN^C?' 2 ) m f° ur redshift bins covering z < 1.4. The 
zero-point of the x-axis has been renormalized to the av- 
erage Q24 of the IR templates at the centre of the redshift 
slice. A value of 2 (0.5) on the y-axis implies that at a 
given value of (724 the relative abundance of AGN and 
SF systems is 2:1 (1:2). Within the errors Iag/n ((7, z) is 
consistent with being unity across the whole width of the 
star forming locus at all redshifts. Deviations from the 
generally smooth variations of i^Q^ l (q, z) can occur on 
the edge of the assessed range of 924 due to fluctuations 
caused by poor statistics. There is weak evidence for 
a gradual decrease of the AGN fraction from about 2 /3 
to roughly 1/3 as one goes from the region which hosts 
sources with radio-excess to that populated by sources 
with excess IR emission. This trend is barely significant 
but interestingly enough it tilts in the opposite direc- 
tion as would be expected if, e.g., AGN activity were to 
manifest itself by exciting increased hot dust emission in 
the MIR. (Note that in general the radio emission could 
also be altered by the presence of an AGN, thus making 
the observed slope less easily interpretable. However, the 
fact that the distribution of f™^ 1 - which has a radio 
contribution that is identical to that in f A Q^ - is essen- 



tially flat, suggests that the radio emission is not strongly 
affected by the AGN.) 

The calculation of f^Q^ n (g, z) involved slightly wider 
probability bins of width l /s (to ensure convergence of 
the distribution function) and was limited to z < 1.1 
due to the ubiquity of 70 /xm non-detections at higher 
redshift (cf. upper panel of Fig. [21]). iJ^Q^iq, z) ap- 
pears to be a constant function of (724, obs with no traces 
of being tilted as detected with marginal significance for 
^AGiT^' z )-> exce Pt maybe in the redshift bin z g [0.72, 
1.06]. Overall, we can thus deduce that our optically 
selected AGN and SFGs occupy the SF locus in very 
similar proportions. A possible explanation for this is 
that both the IR and radio emission are predominantly 
powered by star formation rather than AGN activity. It 
is also conceivable, however, that other (combinations 
of) astrophysical processes conspire to place AGN hosts 
close to the IR-radio relat ion (e.g. I Sanders et al.l 119891 : 
IColina fc Perez-Qledll995l) . 

6.4. Variations of IR/ Radio Ratios with Luminosity 

In a recent work on local IR galaxies fRieke et"aLl (]2009[ ) 
have found evidence of variations in the X-corrected av- 
erage 24 /im/1.4 GHz flux ratio (724.0 with IR luminosity 
Ltir- According to their analysis (724,0 is a constant 
function of luminosity at Ltir < 10 11 L© and then be- 
gins to rise with increasing luminosity. Using the NVSS- 
and IRAS-detected SDSS galaxies, Moric et al. (in prep.) 
see an opposite trend of decreasing FIR/radio ratio when 
they examine (7fir vs. £1.4 ghz for various types of active 
galaxies (both star forming and AGN-bearing). 
We investigate whether or not the star forming sources in 
our sample show any evidence of variations of giro with 
IR or radio luminosity. Since our iiT-corrected monochro- 
matic IR/radio ratios are based on the best-fitting TIR 
SEDs, all luminosity-dependent trends they display will 
be qualitatively identical to those measured for gxiR- 
Comparisons with previous studies are therefore possi- 
ble even if these used a different IR/radio parameter. 
Note that the fact that we are plotting (7tir against lumi- 
nosity implies that upper and lower limits cannot always 
be unambiguously placed along the ordinate. An exam- 
ple are the radio-selected sources in the upper panel of 
Fig. [22] of which we merely know that they must lie to 
the lower left of their limits. They are indicated by an 
arrow pointing diagonally downward. The calculation of 
the median qxiR in a given bin of luminosity should cor- 
rect for measurements that in truth belong to a fainter 
luminosity range. To account for this we construct broad 
luminosity bins (Alog(L) « 1) and assume that most of 
the censored measurements would come to lie in the next 
lower luminosity bin (anything fainter would imply that 
they are more than 3 a outliers to the IR-radio relation). 
We can then 'average' away the effect of falsely assigned 
measurements by (1) computing the median ((/tir) in two 
sets of luminosity bins which are offset by half a bin 
width and then (ii) averaging the two estimates of the 
median thus obtained and reporting the new value half 
way between the centres of the two involved bins along 
the luminosity axis. The medians themselves are calcu- 
lated by applying survival analysis to the jointly IR- and 
radio-selected data as previously done in ^ l6.H[6~3l 
The results of this procedure are shown in the larger 
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two windows in Fig. [2H Using the COSMOS data we 
see no evidence of an increase i n the IR/ra dio ratio at 
Ltir ~ 10 11 L© as suggested bv lRieke et all (|2009l ). We 
do detect a higher value of (qtir) in the brightest IR lu- 
minosity bin but this increase ha ppens around Ltir. ~ 
10 13 L Q , similar to the results of lYounger et al.l (2009). 
It shoul d be mentioned, ho wever, that the methodology 
used bv iRieke et al.l <j2009f ) to derive (724,0 differs signifi- 
cantly from the one used here in that it involves - for 
example - luminosity-dependent (and template-based) 
conversions of IRAS 25 (im flux densities to 24 /im MIPS 
equivalent values. 

While no universal trend for variations of q with IR lumi- 
nosity are detected in our sample we do find that <?tir is a 
decreasing function of radio luminosity (see lower-most 
window in Fig. I22[) . The trend is consistent and in- 
creases rapidly at Li.4GHz ~ 10 24 W/Hz. This could po- 
tentially be the effect of contaminating AGN at high ra- 
dio luminosities in our optically selected sample of SFGs. 
However, the fact that Moric et al. (in prep.) see a 
similar trend in local SF, composite and AGN-bearing 
systems which have been classi fied based on the stan- 
dard optical line em ission ratios (|Kauffmann et al]l2003b 
IKewlev et al.ll2006h suggests that the trend is genuine. 
The two narrower windows in Fig. !22l show the variations 
of the dispersion of Qtir with IR (top) and r adio lumi- 
nosity (bott om). In the low - redsh ift samples of lYun et al.l 
pOOll) and IBressan et all (|2002T) an increase in scatter 
with infrared luminosity is detected. In the present data 
a similar - albeit very weak - tendency is seen; the re- 
duced accuracy of the Ltir measurements of the high-z 
galaxies likely masks most of the trend if present. The 
plot of <7g TIR vs. Li.4GHz, on the other hand shows a clear 
increase in the scatter which starts to manifest itself at 
the same radio luminosity at which the strong decline of 
(^tir) sets in. 

6.5. The IR-Radio Relation at z > 2.5 

While in the previous sections we usually tacitly plot- 
ted data points from high-z sources the fitting of evolu- 
tionary trends in §§ 16. II and 16.21 was restricted to galaxies 
at z < 1.4. This corresponds to the redshift at which 
the 4000 A bre ak leaves the reddest Subaru band with 
deep coverage (|Taniguchi et al.l [20071 ). the z-band. Af- 
ter z ~ 1.4 the break is constrained by the NIR data of 
the J, H and K s bands (McCracken et al. 2009, subm.; 
Capak et al., in prep.). These exposures of the COS- 
MOS field, however, are two magnitudes shallower and 
have gaps between filters, leading to large uncertainties 
in the photometric redshift estimates. Beginning from 
about z ~ 2.5 the Lya (1215 A) break enters the wave- 
length range cover e d by the ground base d photometry 
(jCapak et alj|2007t iTaniguchi et al.l 120071) . As a conse- 
quence the accuracy of the photometric redshift improves 
to again a( Az /(i+z)) ~ 0.03. 

In an assessment of ongoing spectroscopic follow-up ob- 
servations of high-z sources in the COSMOS field, Capak 
et al. (in prep.) find that photometric redshift estimates 
of genuine high-z sources may be scattered to low red- 
shift due to confusion between the Lya and 4000 A break. 
Most of the confusion is due to regions of the Lya forest 
which are not as opaque as expected and/or light from 
nearby foreground galaxies contaminating the apertures. 



Conversely, there is little evidence for any upward scat- 
tering of galaxies at low and intermediate redshift to z > 
2.5. This implies that sources with photometric redshift 
estimates >2.5 represent, with high likelihood, a clean - 
albeit not complete - sample of high-z objects. 
Our sample contains more than ~140 sources at redshift 
z > 2.5, of which approx. 60% have direct detections at 
1.4 GHz and in at least one MIPS filter. As far as we 
are aware, this is the largest sample of high-z sources so 
far, for which it is possible to study the IR-radio rela- 
tion based on direct detections rather than flux limits. 
We must point out, however, that only 2% of the high-z 
sources have a direct detection at 24 and 70 /im while 
the SEDs of the remaining 98% are only constrained by 
a direct detection at 24 fim and an upper flux limit at 
70 fim. Accordingly, the calculated values of Ltir lumi- 
nosities must be regarded as fairly rough estimates of the 
true IR luminosity of these sources as they are primarily 
based on measurements made at a rest frame wavelength 
of -6 /im. iMurphv et~aH (|2009bD caution that the IR lu- 
minosities of high luminosity and high redshift sources 
(Ltir/lq > 10 125 ; z > 1.4) are generally overestimated 
by a factor of ~4 even after subtraction of a flux con- 
tribution from AGN. However, in view of the COSMOS 
study of Kartaltepe et al. (subm.) - who, in the same 
range of IR luminosities, do not see this trend and in- 
stead report that IR luminosities based solely on 24 /im 
data tend to be underestimated in general - we refrained 
from applying any corrections to our data. 
Bearing in mind these uncertainties we plot the 
TIR/radio ratios of our high-z sources in Fig. [23] (left 
panel). For illustrative purposes the measurements of 
Qtir are coloured according to their probability Pr (SF). 
We caution, however, that this classification is based on 
the fiducial (u — K) cut used throughout the paper so far 
and that the evidence presented in Fig. [S] indicates that 
this threshold is no longer appropriate at z > 3. In view 
of this we do not distinguish between star forming sys- 
tems and AGN for the high-z sources but use this global 
sample to derive the average IR/radio flux ratio. The 
right hand side shows the distribution function of qtir 
which is broad (a = 0.59±0.05) and has a median of 
2.71 

-0A4- This value i s in g ood agreement with the local 
measurement of iBelll (|2003l ) (dashed line) and is almost 

identical to the average value of 2.76^oo8 we nn< ^ f° r the 
COSMOS data in our lowest redshift bin in §0 The 
average IR/radio properties of our high redshift sample - 
the most distant sources of which are detected when the 
universe was only ~1.5 Gyrs old - are thus very similar 
to those observed in the local universe. It is important 
to remember, however, that at z > 2.5 the COSMOS 
data contains mostly extremely IR-luminous HyLIRGS 
(Ltir > 10 13 Lq) which are a very different kind of ob- 
ject than those encountered at z < 0.5 where the major- 
ity of our sources have 10 11 < l tir/l q < 10 12 (cf. Fig. 
ED- 

7. DISCUSSION 

Various parameterizations of the IR-radio relation ex- 
ist. The flux ratios ^tir/fir and <Z70 predominantly 
reflect the IR and radio emission of the ISM which is 
caused by two stages in the life cycle of massive stars; 
(i) the main sequence phase during which UV light is 
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converted into FIR emission by dust grains and (it) su- 
pernovae explosions inducing synchrotron emission when 
their shock waves accelerate cosmic ray electrons in the 
galactic magnetic field. The parameter 524, on the other 
hand, is more sensitive to hot dust emissio n triggered by 
AGN activity. Sever al recent papers (e . g. lGarretti | 2002| : 
Appleton et al.ll2004t iFraver et al.ll2006t llbar et alJbooa 
Younger et alJl200ftlMurphv et al.ll2Q09bf ) using 1.4 GHz 
data provide consistent evidence that the local IR-radio 
relation holds out to high redshift. An identical con- 
clusion has been reached using radio flux d ensity mea- 
surem ents at 610 MHz rather than 1.4 GHz (jGarn et al.l 
120091) . 

A majority of the samples previously used to study the 
IR-radio relation in distant galaxies contain only several 
dozen to roughly a hundred galaxies at z < 0.5. The 
COSMOS sample used in the present analysis increases 
the number of sources at redshifts z > 0.5 by at least a 
factor of 5. Also, it is probably the first data set in which 
the number of directly measured IR/radio ratios at high 
redshift is larger than the number of censored values. 
We have studied the IR/radio properties of both SFGs 
and AGN. Various complementary indicators of AGN ac- 
tivity are present in the COSMOS data base. Rather 
than identifying AGN based on a combination of mul- 
tiple parameters we have chosen to work with a single 
statist ical criterion based on the work by ISmolcic et al.l 
( 2008) which classifies radio- and IR-detected sources as 
SF or AGN based on their rest frame optical colour. This 
statistical approach allows us to treat our sources in a 
probabilistic way which is especially advantageous for 
the assessment of systematics inherent in the selection 
of SFGs and AGN. 

7.1. The IR-Radio Relation at Intermediate and High 

Redshift 

At z < 1.4 where both photometric and spectro- 
scopic redshifts have a high accuracy we found no 
compelling evidence of an evolving IR-radio relation. 
The mean IR/radio flux ratios in this redshift range 
are 924,o = 1-26 ± 0.13, gy^o = 2.23 ± 0.13 and 
9tir = 2.57 ± 0.13, where the first two are computed 
using if-corrected IR and radio flux densities. With 
the exception of the highest redshift bin, the median 
IR/radio ratios in the different redshift bins covering the 
range z £ [0, 1.4[ are offset from the averages q by at 
most half the dispersion in the local IR-radio relation. 
An alternative to testing the constancy of rest frame 
IR/radio flux ratios derived using IR SEDs of local 
SFGs, is to form IR/radio ratios with observed flux 
densities (cf. 16.1. 1|) . Variations with redshift of (724, obs 
and g7o, obs in the SF population can then be compared 
with the changes that would be expected for redshifted 
local galaxies (e.g. llbar et al.l[200l iSeymour et"aT1l2009f ) 
in order to assess if the IR/radio properties of the 
latter are compatible with those of high-z systems. Our 
analysis has shown that the decline with redshift of 
observed and predicted IR/radio ratios are indeed in 
good agreement. 

Since measurements of (724 and (770 have been carried 
out with flux-limited data sets which likely are subject 
to some selection band-related bias (see §§ 15.11 and 17. 2[) 
only our value of <?tir lends itself to a stra i ghtfor ward 
comparison with local measurements. iBelll (2003) find 



qriB = 2.64 ± 0.02 for a sample of local SFGs and 
show that this figure is in excellent agreement with the 
IRAS-based FIR measurement ^fir = 2.34 ± 0.01 of 
lYun et al.1 (|2001| ). The COSMOS measurements of the 
mean IR/radio properties of high redshift galaxies thus 
are fully consistent with the local average. Furthermore, 
a subsample of highly redshifted (z > 2.5) COSMOS 
galaxies has a median 2. 71 io 14 which also agrees well 
with both the COSMOS data at low redshift and 
independent local measurements. 

Due to Malmquist bias we detect systems with very dif- 
fering star formation rates (SFR) in the low and high 
redshift universe where only extreme starbursts of the 
HyLIRG class (SFR > lO^gyr" 1 ) are visible. That 
starburst galaxies which we observe when the universe 
was just 10% of its current age follow the same IR-radio 
correlation as local galaxies runs counter to expectation 
which would pre dict higher-than-averag e IR/radio ratios 
for such systems. [Liscnfcld et al.l (|1996l ) have shown that 
a strong (and prompt) enhancement of the magnetic field 
strength is required lest a deficit of synchrotron emission 
develop due to high inverse Compton losses of the cosmic 
ray electrons in the strong radiation field generated by 
the starburst. 

On the other hand, it could be that precisely this inabil- 
ity to detect sources with SFRs of 'normal' (disk) galaxies 
at higher redshift is hiding changes in the IR/radio ra- 
tios of the SF population. In particular, such changes are 
expected as soon as the energy density of the cosmic mi- 
crowave background (CMB) exceeds that of the galactic 
magnetic fields, whereupon inverse Compton losses off 
CMB photons begin to dominate synchrotron emission. 
Based on the typ ical magnetic field strengths in spiral 
arms (a few /xG), iCarilli et al.1 (|2008h estimate that this 
could be the case for regular disks as early as z ~ 0.5. 
Regardless of the remaining uncertainties surrounding 
the maintenance of the local IR-radio relation out to 
high redshift the observational fact per se is an impor- 
tant confirmation of the central assumption used in stud- 
ies that have computed the SF history of the universe 
(lHaarsma et alJl2000HSevmour et aTll2008t iDunne et al.l 
2009: ISmolcic et al.l 12009) using deep radio surveys, 
namely that of an unchanging proportionality between 
SFR and radio luminosity. The constancy of the IR- 
radio relation not only implies that IR and radio mea- 
surements are equally good tracers of star formation out 
to high redshift. It also suggests that the physical pro- 
cesses of massive star formation when the universe was 
only 1-2 Gyrs old used to be strikingly similar to those 
at play in the local universe. 

Currently available observations are not yet capable of 
revealing the exact workings of the astrophysical mecha- 
nisms that cause galaxies at intermediate and high red- 
shift to lie on the IR-radio relation. The improvements 
in sensitivity and spatial resolution with Herschel and 
the extended Very Large Array (EVLA) will be a crucial 
step forward in this respect. A better sampling of the 
FIR emission from distant (U)LIRGs will reveal if their 
SEDs indeed match those of their local analogues as this 
study and previous work generally assume. Similarly, a 
systematic study of the range of radio spectral indices 
in SFGs is important to quantify the inaccuracies that 
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are introduced by assuming a single generic value. Clar- 
ifying both these issues is essent ial if, e.g., trends for a 
luminosity dependence of q (see iRieke etahl 120091 and 
also our Fig. l22j) are to be put on a solid basis. 
Apart from SFGs, the IR-radio relation is also ob- 
served by many AGN-bearing and/or composite 
systems (e.g ISanders et al.l [l98a iMarx et all Il994t 
Sopp fc Alexanderlll99lURov et al"lll998l; iMurphv et al.l 



2009b; ISevmour et al.l l2009h. a lb eit with a larger 
dispersion (iCondon et al.l ~ ll982l IQbric et al.l [2006; 
IMauch fc S adler 20071). In the COSMOS sample we also 
observe that AGN and SFGs often have very similar 
values of (724 and q^. This is true not only for optically- 
selected AGN but also for X-ray detected sources, 
implying that the phenomenon is not merely an artifact 
of our statistical colour criterion for the separation 
between AGN and SFGs. It should be emphasized that 
the finding is not an artificial consequence of template 
fitting; the abundance fraction of AGN and SFGs has 
been derived using observed IR/radio ratios and hence 
does not involve any assumptions about the value of the 
radio spectral index or the shape of the IR SED. 
We find that our optically-selected AGN and SFGs 
populate the locus of the correlation in nearly equal 
proportions out to at least z ~ 1. It is important to bear 
this in mind when IR/radio ratios are used - possibly 
in combination with other indicato rs - to distinguish 
between AGN an d SF sources (e . g. iDonlev et al.1 120051 : 
iPark et al.1 120081 ISevmour et all 120081 ): radio-excess 
outliers can only be used to single out radio-loud AGN 
rather than a complete AGN sample and, conversely, 
selecting only objects that follow the correlation will 
result in a mixture of radio-quiet AGN and SFGs rather 
than a pure SF sample. The fact that an AGN is present 
does not exclude coeval sta r formation in the host galaxy 
(e.g. iSilverman et al.ll2009T ). However, whether or not it 
is responsible for the similar in IR/radio properties of 
SFGs and AGN is still debated, with both supporting 
and contesting ev idence being advanced (e.g. in the case 
of PG QSOs: see ISanders eTaD[l989l: lBarthelll2006l ). 

7.2. Biases Revisited 

Even though the current data set has been selected 
both in the IR and radio, the fact that the average 
IR/radio ratios of the jointly selected sample are usually 
close to those of the IR-selected sample suggests that the 
average values (q) we quote are not the 'intrinsic' value 
one would hope to find in an entirely unbiased sample. 
By separately studying an IR- and radio-selected sam- 
ple of SFGs, however, it is at least possible to bracket 
the unbiased average IR/radio properties. As expected 
the jointly selected sample always lies within this region 
(with the exception of one case - see Table [6]- where the 
amount by which the median of the jointly selected sam- 
ple exceeds that of the IR-selected sample is still smaller 
than the width of the bins used for the calculation of the 
distribution functions of q in the doubly censored data 
set). 

The shift between the average IR/radio properties of 
an IR- and a radio-selected sample are in principle pre- 
dictable based on the dispersion of the IR-radio relation 
and the slope of the differential source counts (cf. equa- 
tion ([5])). While comparable to the locally measured dis- 
persion at intermediate redshift, the scatter of the re- 



lation in our COSMOS sample is significant at z > 1 
where it is twice as large as at low redshift. Part of the 
increase could be due to an intrinsically higher scatter 
a„ at high luminosities (jYun et al.l 120011 : iBressan et al.1 
2002) but probably is also a consequence of the limita- 
tions that are inherent in the calculation of IR luminosi- 
ties when the IR SED is constrained by few points 25 . Un- 
der the simplified assumption of Euclidian source counts 
(/3 = 2.5) equation ([5]) predicts the offset Agbias ~ 0.35 
dex, found in our lowest redshift bins. It also makes a 
fair prediction of a shift of ~0.7 dex between the IR- 
and radio-selected sample at z ^ 1 if one accounts for 
the larger scatter and the finding that at faint fluxes 
IR dChary etail [200l iPapovich et all 12001 and radio 



(|Richarddl2000l; iFomalont et al.ll2006l: iBondi et al.ll2008D 

source counts are sub-Euclidean 26 . 

Figs. [13l [15] and [17] show that apart from biasing the 
average IR/radio ratio, selection effects can also produce 
spurious evolution. Based on the radio-selected sample 
alone we would infer a decrease of the mean (q) out to 
z ~ 1 (see also numeric al values in T ables [3l [5] and [6} . 
As the only recent study, ISevmour et al.1 (|2009h ineasure 
an average value of q-jo for their radio-selected sample 
that is reduced by 0.25 dex at z ~ 1 with respect to 
low redshift. This is close to the evolution d {no)/dz rj 
0.2-0.3 dex which we see for the radio-selected COSMOS 
sample. Based on the argument outlined above a simple 
ex planation for this cou ld be selection effects. However, 
as ISevmour et al] ([2009) adopt the different approach of 
stacking radio sources that are not detected at 70 /im 
rather than including them in the analysis in the form 
of flux limits as we have done, other explanations cannot 
be ruled out. 

As stated in ij l5.ll there is ample evidence from the re- 
sults summarized in Table that the offset predicted by 
equation ([5]) not only occurs in our data but also can rec- 
oncile most apparently discrepant measurements of mean 
IR/radio ratios in the literature. The one exception to 
this generally encouraging agreeme nt are the hi g hly in - 
consistent rad i o stac king results of iBovle et al.l (1 2007) , 
iBeswick et all (|2008l ) and lGarn fc Alexander! (|2009f ) who 
have all studied the mea n (724 as a function of IR flux. 
iGarn fc Alexander! (|2009[ ) in particular pointed out that 
the field-to-field variation of the mean IR/radio ratio can 
be considerable. The prospects are good that the is- 
sue will soon be resolved with the aid of significantly 
deeper EVLA observations at the ~/iJy level that will 
even directly detect the radio emission of the faintest 
24 /im sources. 

8. SUMMARY AND CONCLUDING REMARKS 

On the preceding pages we have discussed the IR/radio 
properties of star forming galaxies (SFGs) and active 

25 Note, however, that in this respect the current data still rep- 
resents an improvement over previous survey-scale samples. 

26 The larger measured scatter cr q tends to increase A%ias while 
the sub-Euclidean counts (with f) ~ 1.5; see references in text) have 
the opposite effect (with respect to a Euclidean slope). Consider 
the different measurements of (qtir) at z £ [1.14, 1.51[ (see third 
row from bottom in TablelrTt as an illustration that the combination 
of the two factors leads to a prediction of Ag^ias which agrees 
excellently with the data. Using equation J5J with /3 = 1.5 and 
u q = 0.78 one finds Ag^ias = 0.70. The observed shift between the 
median (qtir) m the IR- an d radio-selected sample, on the other 
hand, is ~0.67. 
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galactic nuclei (AGN) in the redshift range < z < 
5. Our analysis has benefitted from the extensive multi- 
wavelength coverage of the COSMOS field: each of our 
sources has a flux constraint at 24 /zm, 70 /zm and 1.4 GHz 
from Spitzer/MIPS and VLA observations, multiwave- 
length photometry in ^30 bands from the UV to the 
NIR, and in some cases X-ray detections. Roughly one 
in three of our sources has a spectroscopic, the rest an 
accurate photometric redshift measurement. 
Our primary focus was the evolution of the IR-radio re- 
lation with cosmic time. With a total of approx. 4500 
sources our sample is the largest one which has so far 
been used to study the IR-radio relation at intermediate 
and high redshift. This is particularly true at redshift 
z > 2.5 where we have detected nearly 150 sources of 
which, for the first time, more objects have direct IR 
and radio detections rather than upper flux limits in one 
of the two bands. To our knowledge this is also the first 
time the evolution of both monochromatic and integrated 
IR/radio ratios has been consistently derived using ob- 
jects from the same field. The computation of IR lumi- 
nosities using complete SED template libraries represents 
a further improvement over previous work which has of- 
ten relied on single starburst templates, e.g. that of M82. 
The average IR/radio ratio is subject to selection biases 
in that it depends on the band in which a population is 
selected. We have shown that the average IR/radio prop- 
erties of IR- and radio-selected samples of galaxies differ 
by an amount which is in agreement with theoretical ex- 
pectations and that studies relying on objects selected 
in only one of the two bands run the risk of inferring 
spurious evolutionary trends. For this reason, we base 
our analysis on a sample jointly selected at IR and radio 
wavelengths in order to eliminate biases as best possible. 
Furthermore, we make frequent use of the methods of 
survival analysis which permit us to include all informa- 
tion carried by flux limits from sources that have failed 
to be detected at either IR or radio wavelengths. 
To summarize our findings: 

(i) The median IR/radio ratios of SFGs show little 
variation in the redshift range z < 1.4 and the IR- 
radio relation remains similarly tight out to z ~ 1 
as it is in the local universe (see §§ 16.11 and 16. 2[) . 
Above z ~ 1 the dispersion in the COSMOS popu- 
lation increases, probably due to the reduced accu- 
racy of our measurements and/or an intrinsically 
larger dispersion of IR/radio ratios at high lumi- 
nosities. 

(ii) For a sample of sources with high confidence red- 
shift estimates in the range z > 2.5, the average 
IR/radio ratio is still the same as that found in the 
local universe (see § 16. 5|) . 

(iii) At both 24 and 70 /im many of our optically- 
selected AGN have similar IR/radio ratios as SFGs 
(see § 16. 3|) . The relative abundance of AGN and 
SFGs in our sample is about 1:1. The same applies 
to X-ray detected AGN, of which a large fraction 
has IR/radio ratios which lie in the range measured 
for star forming systems. 

(iv) The median IR/radio ratios of SFGs consistently 
decrease as a function of radio luminosity. On the 



other hand, they remain constant over three orders 
of magnitude in IR luminosity. Only at the highest 
IR luminosities (Ltir ^ 10 13 ^o) nas a tendency 
for an increase been detected (see § I6.4|) . 

(v) Apparently discrepant measurements of the aver- 
age IR/radio ratio reported in the literature can 
be reconciled if one properly accounts for the selec- 
tion band of the respective samples (see § 17.21) . 

In this work we have applied a statistically more sophis- 
ticated treatment of IR/radio ratios than has previously 
been adopted. It accounts for both detection limits and 
selection biases inherent in the data, and our results pro- 
vide firm support for previous reports that the IR-radio 
relation remains unchanged out to high redshift (z > 
4) . The observed constancy is striking evidence that the 
interplay between the life cycle of massive stars and the 
interstellar matter (ISM) has followed a very similar pat- 
tern for more than 10 billion years. IR and radio mea- 
surements apparently represent equally good tracers of 
star formation over much of the history of the universe. 
This has been a critical assumption underlying measure- 
ments of the cosmic star formation history with deep 
radio surveys, and fundamental to the estimation of red- 
shifts for optically undetected sub-mm galaxies. 
The upcoming generation of IR and radio observato- 
ries like Herschel and the Extended Very Large Array 
(EVLA) are expected to be able to perform measure- 
ments of sufficient sensitivity and spatial resolution to 
provide clues on the mechanism shaping the IR-radio re- 
lation in distant star forming galaxies with moderate star 
forming rates. This will be a major advance over the 
mere observation of the phenomenon as is presently the 
case. 
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APPENDIX 
DATA DESCRIPTION 

In the following, additional information on our data sets is provided. Since the determination of source fluxes and 
positions is central to the accuracy of our IR/radio ratios and the band- merging between the observations at different 
wavelengths, we describe the derivation of these quantities in particular detail. 

VLA-COSMOS Radio Observations 

Following the completion of the VLA-COSMOS Large Project ([Schinnerer et al.ll2007f) which had mapped the entire 
2 deg 2 COSMOS field at 1.4 GHz with the VLA in antenna configurations A and C for 275 hours, the central 0.84 
square degrees were re-imaged during an additional 62 hours in configuration A in the spring of 2006. The resulting 
VLA-COSMOS 'Deep Project' mosaic has a resolution (FWHM of synthesized beam) of 2.5" x 2.5" and a pixel scale 
of 0.35"/px. The mean rms sensitivity is <0.01 mJy/beam at the centre of the field and more than 60% (80%) of the 
field has an rms level better than 0.02 (0.03) mJy (see Fig. 3 in Sc hinnerer et al. 2009, subm.). The data reduction 
was carried out with AIPS (Astrono mical Image Processing System; [Greisen 2003) and followed standard procedures 
described in iSchinnerer et alj ((20071 ) and (2009, subm.). 

The AIPS source/component finding task SAD was used to detect radio sources in the Deep Project mosaic in multiple 
iterations with cut-off levels of successiyely lo wer S/N. The resulting source list was then combined with the VLA- 
Large Project catalog ([Schinnerer et al.l 12007) to create a list of 1.4 GHz sources (henceforth referred to as the 'Joint 
catalog') that are at least of 5 a above the local background. Of the ~2900 sources in the Joint catalog 51% are 
found to be unresolved. In these cases the integrated flux density is set equal to the peak flux density of the sources 
which was measured with the AIPS task MAXFIT. In order to correct for bandwidth smearing a position-dependent 
correction based on a model sensitivity map of the Deep Project mosaic has been applied to th e peak flux values. 
Integrated flux densities (which are not affected by bandwidth smearing; e.g. iBondi et al.l [2008) for resolved radio 
sources were determined by integrating over the size of the best-fitting elliptical Gaussian component with the task 
JMFIT. For another 131 (4.6%) of the sources in the catalog which were best fit by the sum of more than one Gaussian 
flux component, the task TVSTAT was used to measure the flux within a manually defined contour around the source. 
The distribution of 1.4 GHz flux density values in the Joint catalog is shown in Fig. 15 of Schinnerer et al. (2009) 
where a detailed description of the construction of the catalog is given. The flux measurements carried out on the 
VLA-COSMOS radio maps have b een compared to those of about 300 so urces also observed at 1.4 GHz in the context 
of the NVSS and FIRST surveys (|Condon et al.lll998t IWhite et al.lll997t ). The agreement between the different data 
sets is reasonable except for a number of NVSS sources where th e VLA-COSMOS observations have probably resolved 
out a large extended flux component fcf. ISchinnerer et al.ll2007T ) . 

For the multiple component sources which often have a complicated radio morphology due to outflows or continuum 
emission from star form ation activity we set the most likely source position by visual inspection of HST ACS images 
(|Koekemoer et aL| [2007) . In all other cases the source position is located at the peak of the radio emission as determined 
by MAXFIT. This definition was also adopted for a small fraction of multiple component sources that cou ld not be reliably 
associ ated with an optical source. The typical accuracy of the radio positions is 0.13" as shown bv ISchinnerer et al.1 

(pool . 

While the Joint catalog serves as the ba sis of our radio-selected sample (introduced in ij |2.3.1[) we will also use the 
MIPS 24 fim catalog described in § IA.2.1l to construct a sample of infrared-selected sources. In addition to counterparts 
with S/N > 5 from the Joint catalog we also allow for sources having S/N > 3 in the presence of IR detections. To 
find additional sources in this S/N range we ran MAXFIT at the position of IR detections without a counterpart in 
the Joint catalog. The MAXFIT box size is chosen in accordance with the uncertainty in the position of the IR source 
(see Appendix [Bj). Since the vast majority of sources at these low d etection level s are l ikely unresolved (radio sources 
with flux density <0.1mJy typically have sub-arcsecond sizes, see IBondi et al.l ([20081 ) and references therein) their 
integrated flux is equated to the peak flux density. 

If the extracted peak flux density has a significance level of less than 3 a or if MAXFIT fails to converge, upper limits 
on the total flux from the local sky background are used as flux constraints for a potential source at that position. 
The noise in the Deep Project image was estimated with the AIPS task RMSD in a box of dimensions 105" x 105" that 
was moved across the survey area in steps of 2.45" in R.A. and Dec. We adopt an upper flux limit of 3 times the 
local rms noise level for the rest of the analysis for undetected radio sources. The numeric value of this upper bound 
is obtained by reading the noise image at the pertinent position. 

S-COSMOS Infrared Observations 
MIPS 24 [im Data 

During MIPS Cycle-2 and Cycle-3 the S-COSMOS project ([Sanders et al.ll2007h imaged the whole COSMOS field at 
24 /im in medium and slow scan mode, respectively. The data fr om both MIPS cycles wa s coadded and combined with 
the MOPEX package (MOsaicker and Point source Extractor; iMakovoz fc K han 20QE) after elimination of transient 
sources such as asteroid. The resulting mosaic has a pixel scale of 1.2"/px, while the FWHM of the MIPS 24 fim 
PSF is 5.8". About 90% of the area was mapped with a median integration time of ~3400 seconds resulting in a 
la sensitivity of ^0.018 mJy. Over the remaining ^10% of the field the average integration time was roughly 7000 
seconds, leading to an equivalent depth of 1 a ~ 13 ^Jy. A noise map for the S-COSMOS 24 /im observations was 
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generated using the associat ed coverage map. The details of the MIPS 24 /im data reduction and source extraction 
procedures are spelled out in lLeFloc'h et alT (|2009f ) . In the following we summarize the points that are most relevant 
to our work with respect to the 24 /i,m source catalog. 

Source detection in the 24 /im map was carried out with SExtr actor (iBertin fc Arnoutsll 19961 ) which returned positions 
that served as input to the PSF-fitting algorithm DAOPHOT (|Stetsonlll987fl . DAOPHOT performs simultaneous PSF 
fitting to multiple sources which is an important advantage for the deep and crowded MIPS 24 /im images where objects 
may be blended. By inserting and re-extracting artificial sources in exactly the same manner as true astronomical 
sources, iLeFloc'h et al.l (|2009l ) found that the approach of PSF fitting is reliable down to a flux density of ~60 fiJy 
which in most regions of the mosaic corresponds to a S/N of about 3. At this detection threshold there are ~50,000 
sources in the area covered by the 1 .4 GHz observat i ons. Given the resolution of ~6" of the MIPS images most of these 
are not resolved. As described in ILeFloc'h et all (|2009l ). the fidelity of the S-COSMOS 24 /im flux calibration was 
ascertai ned by checking tha t the K s — [24] colours for stars in the COSMOS field listed in the catalog of the 2MASS 
survey (|Jarrett et alJ l2000l) did not deviate from the expected relation. Note also that total flux measurements at 
24 /im account for all the flux in the extended wings of the PSF. The astrometry of the MIPS data is usually accurate 
to within a few tenths of an arcsecond which was confirmed by the cross-correlation with the 2MASS catalog, which 
revealed only a small systematic offset of 0.3" in declination that was subsequently corrected in the catalog of 24 /im 
sources. 

Since the 24 /im catalog extends to a lower detection threshold of ~3 a we do not attempt to extract extra sources 
in the vicinity of radio detections. If no 24 /im counterpart to a radio source is found during the band-merging (see 
Appendix [Bj wc directly assign an upper 3 a flux limit based on the value of the noise map at the according position. 
The conversion from surface brightness noise (as specified by the noise image) and point source detection noise (as 
required for a flux limit) was derived by scaling the median value of all pixels with an exposure time of ~3400 seconds 
in the noise image to the average la sensitivity level of ~0.018mJy. The adopted value of the conversion factor is 
2.3mJy(MJy/sr)~ 1 . 



MIPS 70 nm Data 

MIPS 70 /im o bservations of the COSMOS field were carried out in parallel with the 24 jxm imaging using the scan 
mapping mode (jFraver et al.l 120091 ). The FWHM of the data is 18.6" and the pixel scale 4"/px. As in the case of 
the 24 /im survey a limited area has a high coverage by repeated observations with a total exposure time of ^2800 
seconds. The median effective exposure time over the remaining 90% of the survey area is 1350 seconds leading to an 
average point source noise (la) of 1.7 mJy The 70 //m observations were reduced with the SSC pipeline tools GeRT 
and MOPEX according to standard procedures for MIPS-Ge survey data. Special attention was given to the creation 
of an accurate noise image which represents the uncertainty owing to both small scale (scatter in repeated observations 
of each sky pixel) and large spatial scale noise propertie s (i.e. pixel-to-pixel correlated noise and confusion noise) . All 

steps leading to the final data products are described in lFraver et alJ (|2009l ). 

Source detection and extraction was performed within the MOPEX package (jMakovoz &: Marleaul 120051 ) using the 
APEX (Astronomical Point source EXtraction) peak algorithm, augmented by additional specialized scripts. The 
noise image was employed to detect and then fit peaks with S/N > 3 using the point source response function (PRF) . 
The final source list has been cleaned of spurious detections along the first Airy ring of the bright 70 /im sources. We 
checked all objects down to the flux level at which it would take a >2.5er noise fluctuation in order for a spurious 
source to be flux-boosted to S/N > 3 in the first Airy ring (the amplitude of which is about 2% of the peak). 
The final list of 70 fj,m sources contains almost 3000 entries in the region of the S-COSMOS survey that overlaps with 
the VLA-COSMOS area. Flux measurements in the catalog have been correc ted for the additional flux outside the 
PRF image (i.e. beyond the first Airy ring) using Spitzer Tiny Tim models (jKristl [20021 ) and placed on a constant 
vS v scale. This colour-correction is accurate within 2% for a wide range of galaxy SEDs across the filt er bandpass. 
The calibration of the 70 /im data was confirmed to be consistent with the official MIPS calibration (jFraver et alJ 
2009) which is accurate to 7% (| Gordon et al.ll2007l ). Finally, a small positional offset in declination was found in the 
comparison of radio and 70 /im source positions. However, at ~0.2" it is significantly smaller than the scatter measured 
for the positions of individual sources which amounts to about 2". 

To provide upper 3 a flux limits at 70 /im whenever the band-merging with the radio or 24 /im fails to identify a 
70 /im counterpart we convert the local surface brightness noise estimate given in the noise image to a point source 
noise. After accounting for the flux in the extended wings of the PSF the corresponding conversion factor is 14.9 mJy 
(MJy/sr)" 1 . 



Multi- Wavelength Photometry 

The COSMOS photometry catalog is an i-band selected catalog with PSF-matched photometry (FWHM = 0.6") from 
30 broad, medium and narrow band filters extracted at the positions of the Subaru i + band detections. The wavelength 
range covere d by these observations extends all the way from the UV at 1550 A to the MIR at 8 /im. iCapak et alJ 
(|2007l l2008f) provide a complete description of the observations and data reduc tion leading to the compilation of the 
multi-wavelength data set. Here we use the photometry catalog compiled by lllbert etldl (|2009aD which lists more 
than 600,000 COSMOS galaxies with i + < 26 detected in a region roughly contiguous with the area covered by the 
VLA-COSMOS observations. 
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X-Ray Observations 

The COSMOS field has been obser ved with XMM-Newton for a total of ~1.5 x 10 6 sec onds (—400 hours), resulting in 
a homogeneous depth of -5 x 10 4 sec ([Hasinger et al.ll2007t iCappelluti e"t~al|[2007l |2009) over 1.92 deg 2 . The associated 
catalog includes 1887 point-like sources chosen to have a high probability of being a reliable detection 27 in at least one 
of the soft (0.5-2 keV), hard (2-10 keV) or ultra -hard (5-10 keV) bands d own to limiting fluxes of -5xl0" 16 , ~3xl0 -15 
and — 5xl0~ 15 ergcm~ 2 s _1 , respectively fsee ICappelluti et al.l (120091 ) for additional details). Moreover, the central 
part of the field has been observed with Chandra ([Elvis et al J 120091 ). providing precise positions for about half of the 
XMM sources. Brusa et al. (in prep.) have determined optical counterparts to more than 90% of the XMM catalog 
sources. In th e present work we u se only those with unique/secure optical counterparts and a successful fit to the UV 
to MIR SED ([Salvato et al.ll200"9T ). The optical position has been used to associate our radio and IR sources to XMM 
detections using the same search radii as adopted for the band-merging of our radio and IR sources with the optical 
photometry (see Appendix [5J . 

Redshift Information 

The catalog of lllbert et al.l (|200 9a) contains photometric redshift estimates for all tabulated objects with the ex- 
ception of (i) sources lying in a 'masked' region of the optical imaging due to a bright star in their vicinity, or (ii) 
sources that are detect ed at X-ray wavelen gths, in which case a special set of AGN-templates was employed to derive 

a photometric redshift (Sal vato et al.|[2009l ). 

Spectroscopic data has already been gathered for more t han 20,000 sources in the COSMOS field le.g. lLillv et aHl2007t 
iTrump et al.l 120071 : ITrump et al.ll2009c iLillv et alJl2009b Kartaltepe et al., in prep.; Salvato et al., in prep). Whenever 
possible, we give precedence to redshift information from spectroscopy. In all other cases we assign a photometric 
redshift to the sources in our radio- and IR-selected samples. The details of the band-merging of the radio and IR 
catalogs with the optical data are presented in Appendix [B] while the next two paragraphs summarize the properties 
and screening of the available redshift information itself. 

Spectroscopic Redshifts 

At the time of writing the data base of reduced spectroscopic observations in the COSMOS field included more than 
14,000 sources. Most of these (75%) were observed by VIMOS in the context of the zCOSMOS survey, observations 
with DEIMOS and IMACS(3) account for another ~12% each and a small fraction of sources has LRIS, SDSS or 
FORS1 spectroscopy. In some cases spectra of the same source have been obtained by more than one of these 
instruments. If these measurements do not agree within experimental error we disregard the one with the lower 
quality flag. The quality of the redshift measurement for our spectroscopically observed optical counterparts is judged 
by the spectroscopic confidence classes 28 of the zCOSMOS survey. In particular, we accepted only those objects with 
confidence flag 3 and 4 (regardless of the decimal place, see below), as well as 2.5, 2.4, 8.5, 9.3, 9.5 and 1.5 29 . The 
decimal place indicates the level of agreement between the spectroscopic and photometric redshift measurements. Our 
choice of accepta ble flag values ensures that the reliability of all measured spectroscopic redshifts is better than 99% 
([Lilly et al.H2009h . 

The classification with decimal places was not available for some of the sources observed by follow-up campaigns other 
than zCOS MOS. In such c ases w e upd ated the spectroscopi c confidence flag as required by comparison with the photo-z 
estimates of lllbert et al.l (j2009al ) and iSalvato et a l. (2009). For spectroscopic targets with unsatisfactory confidence 
flags (<3% of the spectroscopic targets) we checked whether a reliable photometric redshift estimate was available 
according to the selection criteria described in the next section and were thus able to recover redshift information for 
all but 4% of these. In Table [1] we summarize the availability of redshift information for our sources. 

Photometric Redshifts 

Sources successfully matched to an entry in the COSMOS photome try catalog generally have a precise photometric 
redshift estimate derived with the code Le Phare dlb ert et al.ll2009al ). The average number of bands available for the 
computation of photometric redshifts for the matched radio and 24 (im sources was 29. The accuracy o f the photometric 
redsh ifts was calibrated with the help of more than 4000 high-confidence redshifts from zCOSMOS ([Lilly et all 120071 
2009) sou rces with i < 22.5 a nd verified using faint 24 /jm sources (Kartaltepe et al., in prep.) in the range i > 22.5. At 
z < 1.25 lllbert et alJ ~ i|2009af ) found that the dispersion a( Az /(i+z)) of photo-z measurements is 0.007, 0.013 and 0.051 
in the Subaru i^ B -band magnitude ranges of < 22.5, 22.5 < < 24 and i~^ B > 24, respectively. Beyond redshift 
z — 1.25 the statistical accuracy of the photometric redshifts abruptly decreases by a factor of —3 with respect to the 
range z < 1.25 because the Balmer break is redshifted into the NIR filters which have less sensitive photometry and 
non-contiguous wavelength coverage. We use these values of the photo-z dispersion to remove sources in our sample 
with uncertain photometric redshifts; the criterion which each source must satisfy to remain in the sample is that its 

27 The maximum likelihood threshold ensures that a catalogued XMM source has a probability of at most ~5xl0~ 5 of being a spurious 
detection. 

28 http : //archive . eso . org/cms/eso-data/data-packages/zcosmos-data-release-dr2/index_html#release_notes 

29 Also Type 1 quasars (i.e. broad line objects), to the flag values of which 10 has been added, and objects serendipitously covered by 
the slit of another spectroscopic target (the confidence flag of such objects are marked by a prepended '2') were admitted into the sample 
as their statistical reliability is equal to the classes listed in the body of the text. 
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photo-z error be smaller than 2 a of the dispersion at a given magnitude and redshift. We also eliminate sources with 
a Xphz above a threshold corresponding to two standard deviations in an ideal ^-distribution (with the number of 
degrees of freedom equal to the difference between the number of filters used for deriving the redshift and the number 
of free parameters, namely three, i.e., the redshift, the template type and the template normalization). In practice 
very few sources (<0.5% in the combined radio- and IR selected samples) are rejected due to the latter criterion. 
A somewhat larger number has been excluded due to a broad redshift probability distribution but they nevertheless 
represent only ~1% and ~2% of the radio- and IR-selected sample, respectively. The vast majority of sources without 
spectroscopically constrained redshift have a photometric redshift estimate within the statistically expected accuracy 
(46.4% and 55.5% in the radio- and IR-sele c ted sam ples, respectively). 

The template library used by lllbert et al.l ()2009af) consists of star-forming and passive galaxy SEDs but includes 
none t hat reflect the features expected in spectra of sources with a dominant contribution from an AGN. iSalvato et al.l 
(2009) have computed photometric redshifts for all XMM-detected sources in the COSMOS photometry catalog using 
templates with varying AGN contributions. Their redshift measurements account for source variability and achieve 
an accuracy of er( Az /(i+z)) < 0.015 at i\ B < 24.5 for both type 2 and type 1 AGN and QSOs out to z ~ 4.5. As no 
confidence intervals were available for the best-fitting AGN photo-z estimates and since the statistical dispersion is 
far smaller than the width of the redshift bins studied below we have kept all sources with redshifts derived from the 
AGN-template library in the sample. 

BAND-MERGING 

Definition of the Reference Position 

• IR-selected sample - The ~6" FWHM of the 24 /im PSF and Spitzer pointing uncertainties result in a non- 
negligible uncertainty on the source centroids of 24 /im sources. In order to have a more precise reference for the 
positional matching with other catalogs we searched for IRAC counterparts to each 24 /jm source (the FWHM 
of the IRAC PSF is about 1/3 of that at 24 /im). Since the IRAC imaging of the COSMOS field is shallower than 
that performed at optical wavelengths but still sufficiently deep to detect a counterpart for most 24 /an sources 
this approach simultaneously reduces the likelihood of assigning the wrong optical counterpart in the subsequent 
band-merging with the optical photometry catalog (see last paragraph of this section). 

We correlated the 24 /jm catalog with the IRAC catalog (|Sanders et alJ2007t lllbert et al. 2009b) using a matching 
radius of 2" and found 45,827 unambiguous counterparts (92% of all sources in the 24 /mi catalog). In 572 cases 
more than one IRAC source was found within the search radius. Unless the closest of the potential counterparts 
was at least twice as close to the 24 /mi position as the other candidates these sources (1% of the catalog) 
were excluded from the sample. In the case of 3668 24 /mi sources for which no IRAC counterpart could be 
determined we searched the photometry catalog for optical counterparts within 0.6", the FWHM of the ground- 
based imaging. If this match was ambiguous the source was disregarded. Unique optical counterparts were 
found for another 1386 sources, however, and the final percentage of 24 /im sources with reference coordinates 
thus increased to 94%. 

• Radio- selected sample - The 1.4 GHz source positions are accurate to within about 0.1", which is sufficiently 
accurate that they can be used directly as the reference position for the identification of optical and IR counter- 
parts. 

Search for 24 l^m Counterparts (Radio- Selected Sample) 

Starting with the radio source position we identified 2233 objects (77% of the radio catalog) that have a single 
24 /im counterpart within 2" while 647 (22%) radio sources remain unmatched. Based on their radio properties the 
24 /im non-detections are primarily tbd . For <1% of the radio-selected sample the match was ambiguous. This 
subset includes two pairs of radio sources which have been assigned the same 24 /im counterpart and 17 for which 
more than one IR-detection was present within the search radius (mostly due to a complicated 24 /tm-morphology or 
the occasional over-deblended source). 

Search for Radio Counterparts (IR-Selected Sample) 

Since the identification of 24 /im counterparts to 1.4 GHz sources in the immediately preceding paragraph was based 
on the precise radio positions these pairs of objects can be taken over directly into the IR-selected sample 30 . Note that 
subsequently we search for additional radio detections with S/Niagbz > 3 in the vicinity of all 24 /im sources without 
a radio counterpart in the radio catalog described in ^ IA.11 The box in which we check for these additional faint radio 
sources has dimensions 1.75"xl.75" (i.e. 5x5 px 2 at the pixel scale of the Deep Project image) which corresponds 
to the FWHM of the PSF of the first two IRAC bands. Radio fluxes for an additional 639 24 /mi sources could be 
measured in this way. 



30 If the according 24 (im source has no trustworthy IRAC or optical counterpart, the 24 /^m-radio association in question is not included 
in the sample. This affects only a very small fraction of the IR-selected sample since the fraction of sources for which no reference coordinate 
could be determined is <1%. 
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Search for 70 fim Counterparts 

We allow for a maximal separation of 6" between reference coordinates and an accepted 70 /zm source. Due to the 
broad FWHM (18.6") of the MIPS 70 fim PSF there is a relatively large probability that the emission from more 
than one 24 fim or 1.4 GHz source is blended into a single 70 /jm source with no distinguishable secondary peaks by 
which the individual components could be separated. This occurred for ^7% of the 1.4 GHz sources with a potential 
counterpart at 70 ^m. In the IR-selected sample 21% of the matched sources had a potential 70 /mi counterpart that 
lay within less than 6" of at least one more 24 ^m source. Unless these ambiguous 70 /im associations are at least 
twice as close to the nearest 24 fim or radio source as to the next best candidate they have been excluded for the 
analysis of the paper. Ultimately, this was the fate of 70% (50%) of the radio (24 fim) sources with blended 70 fim 
photometry. In the interest of achieving a high accuracy of the measured IR/radio properties we prefer this approach 
to, e.g., attempting to fit the 70 fim flux distribution with two components with peak positions fixed to the reference 
position. 

To summarize, in we have found ^1500 reliable 70 fim counterparts to 24 fim sources. A further 60 sources in the 
IR-selected catalog had more than one counterpart within 6" and were removed together with the 24 fim objects that 
are associated with a blended 70 fim detection. The majority of the 24 fim sources (67%) has no directly detected 
counterpart. In the radio-selected sample we were able to identify ^820 unambiguous 70 fim counterparts while 
excluding ^60 sources because they either were associated with a blended 70 fim source or had an ambiguous match. 
For some 70% of the radio sources no 70 fim counterpart could be found. 

Search for Optical Counterparts and Spectra 

A search radius of 0.6" is used for the cross-correlation of the radio-selected sample with the COSMOS photome- 
try catalog and the spectroscopy catalog. This figure corresponds to the FWHM of the PSF-matched ground-based 
photometry. Whenever more than one optical counterpart is found within this distance the according radio source is 
removed from the catalog. This choice was made because any photometric redshift derived for the source would nec- 
essarily be unreliable as both sources contribute to the measured flux. The same applies to spectroscopic observations 
which have a slit width of ~1" and are thus treated analogously because sources within 0.5" could potentially lead to 
ambiguous spectral features. 

To find optical counterparts and spectroscopic information for sources in the IR-selected catalog we use a search radius 
of 1" which accounts for the somewhat larger uncertainty in the predominantly IRAC-based reference positions. Also, 
we ease the rejection criterion for ambiguous matches somewhat in that we accept optical or spectroscopic counterparts 
which are more than twice as close as the other candidate(s), provided the next best candidate lies beyond 0.6". 

DERIVATION OF DISTRIBUTION FUNCTIONS WITH SURVIVAL ANALYSIS 

Doubly Censored Data 

To derive the underlying distribution /(g) of IR/radio ratios q for a sample in which measurements of q can be 
limited both from a bove and below 31 due to non-detections either in the IR or radio band we use the method of 
iSchmitt et al.l (|1993f ). Their approach is based on maximum likelihood estimation and requires no assumptions about 
the form of the true distribution of q (it is thus said to be non-parametric) . If represented in the form of a cumulative 
distribution, F(q), the estimator is equivalent to the probability of measuring a value of q smaller than a given ratio 
(which we will denote by a capital Q): 

F(q) = P(q < Q) . 

In the following we summarize the steps involved in the construction of F(g). 

To begin with, the n data points Qi are arranged in M bins 32 with ascending bin centres Qj > Qj+i, j — 1, M. 
Let dj be the number of detections in the j-th bin which spans the range [Qj — dQ/2, Qj + dQ/2[. For the j-th bin we 
also define Uj , the number of upper limits smaller than Q j and, in analogy, lj as the number of lower limits exceeding 
Qj. Note that in view of the definition of Uj a nd U it is useful t o choose bin boundaries offset by —dQ/2 and +dQ/2 
for the upper and lower limits, respectively fcf. IAvni~et al. 1980]). This eases the computation of the vectors d, u and 
I which can be easily constructed as slightly shifted histograms of the detections as well as of the upper and lower 
limits. 33 

31 The formalism summarized here assumes a random censoring model in which measurements are censored independently of their 
location in the distribution. 

32 Since the bin width can always be decreased to the point that each bin contains only one (untied) detection this assumption is not 
restrictive. 

33 An important modification of the vectors d and u is necessary if the first uncensored value, qx , occurs at a higher value of Qj than 
the smallest upper limit in the data set. By construction, the algorithm of Schmitt ct alj i 19931) is not sensitive to upper limit s th at lie 
below the smallest uncensored data point. As a consequence the normalized cumulative distribution function defined in equation |)C3|) rises 
above for the first time in the bin containing the smallest direct detection. The correct asymptotic value of the distribution function for 
q — > min(Qi), however, should not be in this case but a number larger than or equal to #H, Qi < qx}/ n - Here Qi < qx} is the 
number of upper limits smaller than qx and n is defined immediately following equation ||C2[| . The correct behaviour of the cumulative 
distribution function can be achieved by setting do = Qi < qx} and Uj = for all bins j with bin centres Qj < qx- 

Note that even when there are lower limits that exceed the largest uncensored value, no analogous fixes are needed to correct the doubly 
censored distribution function in the limit q — > max(<2;). 
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The value fj of the distribution function /(g) in the j-th bin is constrained by the probability 

h (CI) 



' 3 



M \ " 3 I M 

3=1 \fe=l / \k=] + l 



of obtaining the triplet (dj, Uj, lj). Intuitively equation (|C1[) can be understood as the product of the a priori 
probabilities of a measurement in the j-th bin being a detection (first term), an upper limit which - if detectable - 
would in truth be loc ated in a lower bin ( second term) or an upper limit in actually located in a higher bin (third 
term). ISchmittl (|1985l ) and lCampbelll (|1981l ) have shown that by introducing the likelihood and by using the constraint 
that X)j=i fj = 1 tlie- probability in equation (jCll) can be used to derive a system of M coupled fixed point equations 
for each of the fj : 

'■^^W^W' <c2) 

k=j + l \<Lim=lJ m / k=l \<Lim=k+lJ m / 

Here n = j (dj + Uj + lj) is the total number of measurements in the sample. Equations (|CJ2[) can be solved by 
iteration starting with fj — dj jn and Uj = lj = and the cumulative distribution constructed according to 

Fj = J2fk, j = l,...,M (C3) 

fc=i 

once convergence has been achieved. This was typically the case after approximately 30-70 iterations. 
Various approaches to calculate th e uncertainty asso ciated with the doubly censored cumulati ve distribution fu nction 
F(q) have been proposed (see, e.g. iZhu fc Sunll2007f ). Here we use the analytical formalism of iTurnbulil (|1974f ) which 
is based on the Fisher information matrix J with elements 34 

t _ dj dj+i lj Uj_ 

" (Fj-Fj.x)^ (F j+1 -Fj)^ (1-Fjf F?' j 1 \1 1 

dj+i 

3,3 + 1 — Jj + l,j — /J7 p« ' 

d_M L Im _ UM_ ( C4 ) 
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J3J+1 - J 
Jmm = 



{F m -Fm-i) 2 (1-Fm) 2 F 2 m 

Jij = 0, for |j-j|>l. 

J is a symmetric Jacobi matrix in which only the diagonal and off-diagonal elements are non-zero. In the limit of 
large samples its inverse, V, is an asymptotically unbiased estimate of the covariance matrix of the Fj (j = 1, M) 
such that one has 

4. = var^) = Vl . (C5) 

Following Fcigclson & Nelson (1985) and references therein it is then possible to define 2 a confidence intervals of the 
form 



F e X p( ± 1.96 CTFj ) ^ (C6) 



which always lie in the permissible range [0, 1]. 



By applying the algorithm of equation (IC2[) to singly cens ored data sets it was possible to check whethe r its output 
agrees with that obtained using the formalism of Appendix IC.2I We found that the v arian ce in equation (|C5[) needed 
to be doubled in order to achieve consistency of the confidence intervals in equation (|C6j) with those computed with 
the algorithm specifically designed to deal with singly censored data. For the sake of consistency the uncertainties 
on the doubly censored dis trib ution functions stated throughout the paper have thus been computed using twice the 
nominal value of equation (|C5|) . 

It can be useful to convey information on the shape of a (cumulative) distribution in form of percentiles. Given Fj 
and the associated confidence intervals, the p-th percentile is defined as that q for which F(q) = p/ioo. Similarly, the 
points where the curves defined by the upper and lower confidence intervals - given by equation (|C6j) - of F (q) are 
equal to p/ioo are used as a measure of the uncertainty on the sought percentile. 

34 If Im = (and hence Fm = 1) the error on Fm is zero by definition. In this case is useful for computational reasons to reduce the 
dimensionality of the problem by one by setting Im—1 = Im— 1 + n M an d thereafter only considering the (M — l)-dimensional vectors d, 
u, I and F. 
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Singly Censored Data 

A formalism which is identical to or simplified with respect to that of the previous section can be used when dealing 
with single censoring. If, for instance, we consider a radio-selected sample its distribution of q is censored from the 
left and can be represented as 



0.1 = max(q i , g*) 

-1 if Qi < q- 



(G7) 



S,= 



if q t e [q t , <?, +1 



in ana logy to equati o ns . The treatment of singly censored data has been discussed in detail infFcigclso n k, Nelson! 
(|1985() and iSchmittl (|1985l ) with an emphasis on astronomical applications. The former paper also describes how to 
transform left censored into right censored data, the treatment of which has been documented more frequently in the 
literature for historical reasons (see references in the two aforementioned articles). Here we simply would like to point 
out that equations (|CJ2|) , which constrain the distribution function of doubly censored data, remain valid in the case of 
single censoring. However, the fact that either all of the Ui or U as defined in § IC.ll equal zero allows the coupled system 
of equations to be solved analytically in terms of the cumulative distribution function Fj throug h the Kaplan-Meier 
produ ct limit estimator (|Kaplan fc Meier! 1 19581) . To do so we follow exactly the prescriptions of Feigclso n fc Nelsor] 

DERIVATION OF THE RELATIVE AGN FREQUENCY F AG n(Q,Z) 
The relative frequency of AGN at a given value of 524 [70] (and in any of the considered redshift slices) is computed as 

, N V - Pr (AGN)- , . 

fAGN ^>- £pr( SF)l ' ™ 

where the summation is carried out over all objects i in the given redshift bin that lie in the pertinent range of IR/radio 
ratios <Z24[7o]- Furthermore, we have by definition that Pr(AGN) = [l-Pr(SF)], as illustrated in Fig. [5J Note that 
if the probabilities of SF- and AGN-'hood' are discretized (e.g. drawn from a set of integer probability values in per 
cent) the numerator and denominator in equation (jDll) could be rewritten as a summation over probabilities. In this 
case the denominator, for instance, takes the form: 

100 

^Pr(SF) ! = ]T A l% xz/100. (D2) 

i i—1 

Here stands for the number of objects with a probability of i per cent of being a SFG. If a data set with only 
uncensored measurements is being studied the computation of fAGN (<7j z) is straightforward since the numbers iVj% 
can be directly determined by counting the suitable objects in each bin of q and z. In the presence of censored data, 
however, iVj% should be regarded as an effective number of objects which accounts for the fact that (i) the true value 
of q associated with a limit found in a specific bin, j, might lie outside that bin (if its actual value were determined 
with the help of better data), or, (ii) limits detected in neighbouring bins have a finite probability of ending up in the 
j-th bin (e.g. lower limits populating a bin with rank smaller than j or upper limits in a bin centred on a larger value 
of q than that of the j-th bin). 

To compute this effective number of objects in the j-th bin, iVj% e a.C?)> we derived the distribution function of objects 
with the desired probability of being SFGs following the techniques presented in Appendix IC.ll Due to the finite 
number of objects in each redshift bin it was necessary to work with a limited number of probability bins to ensure a 
well-behaved estimate of the distribution function fi%{q). Since fi%(q) specifies the fractional contribution of a bin in 
q to the total number of measurements we have that 

N i%< off. (j) = fj, i% x JVtot, i% , (D3) 

where A toti ; % is the total number of objects used to compute the distribution function. With these modifications the 
new expression for fAGN (q, z) in the j-th bin (to be compared with equation (IDip ) becomes: 

f („ A £r=i^,cff. x(l-Pr(SF)) t 
Ei=i A^eff. x (Pr(SF))i 

Note that because of the coarse binning in probability space i no longer defines a number in per cent in this case but 
a range of probabilities instead. Consequently, we also have introduced a new weighting factor, namely the median 
value of the prob abilit ies in the grouped data, (Pr(SF))j. It replaces the probability of the individual objects, i /ioo, 
used in equation (ID2|) . In practice the median usually does not differ much from the mean, i.e. it might take on a 
value of ~0.16 for objects with probabilities in the range 0.1 < Pr (SF) < 0.2. 
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TABLE 1: Summary of sample properties. 



category 


no. of src. 


% of cat. src. 


remarks 


radio-selected sample 


all 


2901 


100 


initial 1.4 GHz sample 


no redshift 


748 


25.8 


includes 21 1.4 GHz sources outside area covered by multi-A photometry 


spec-z 


769 


26.5 




photo-z 


1232+114 


46.4 


regular + AGN photo-zs 


unreliable redshift 


38 


1.3 




useable 


2020 


69.6 


no ambiguous counterparts, reliable redshift & AGN/SF classification available 


Pr(SF) > 0.5 


766 


26.4 


optically selected SFGs with (u — K) < 2.36 (38% of 'useable' sample) 


Pr(SF) < 0.5 


1254 


43.2 


optically selected AGN with (u — K) > 2.36 (62% of 'useable' sample) 


IR-selected sample 


all 


5301 


100 


initial 24 (im sample restricted to S„(24/im)> 0.3 mjy 


no redshift 


1003 


18.9 


includes 304 24 fim sources outside area covered by multi-A photometry 


spec-z 


1254 


23.7 




photo-z 


2712 + 231 


55.5 


regular + AGN photo-zs 


unreliable redshift 


101 


1.9 




useable 


3259 


61.5 


no ambiguous counterparts, reliable redshift & AGN/SF classification available 


Pr(SF) > 0.5 


1822 


34.4 


optically selected SFGs with (u — K) < 2.36 (56% of 'useable' sample) 


Pr(SF) < 0.5 


1437 


27.1 


optically selected AGN with (it - K) > 2.36 (44% of 'useable' sample) 


jointly IR- and radio-selected sample 


all 


6863 


100 


1560/3960 sources selected only at 1.4 GHz/24 /im; 1341 in both catalogs 


no redshift 


1620 


23.6 


includes 321 sources outside area covered by multi-A photometry 


spec-z 


1486 


21.7 




photo-z 


3362 + 275 


53.0 


regular + AGN photo-zs 


unreliable redshift 


120 


1.7 




useable 


4454 


64.9 


no ambiguous counterparts, reliable redshift & AGN/SF classification available 


Pr(SF) > 0.5 


2215 


32.3 


optically selected SFGs with (u-K)< 2.36 (49.7% of 'useable' sample) 


Pr(SF) < 0.5 


2239 


32.6 


optically selected AGN with (it - K) > 2.36 (50.3% of 'useable' sample) 



TABLE 2: 

Representative studies (ordered 



ACCORDING TO THE INVESTIGATED IR/rADIO PARAMETER.) ON THE SPATIALLY UNRESOLVED IR-RADIO CORRELATION USING IR. AND 1.4 GHz DATA. 



survey /field 



reference 



selection band 



other flux limits 



<</} 



HDF-N 
SWIRE 
xFLS/SWIRE 

IRAS 

Subaxu-XMM-New- 
ton Deep Field 

xFLS 

xFLS 

13 H XMM- Newton/ 
Chandra Deep Field 

IRAS 

HDF-N 

var. survey fields 

GOODS-N 

xFLS 

EGS/FIDEL 
misc. literature 



local 

* £ 1 

* £ 2 
z < 1 

2 < 3 

z < 0.16 
z < 1.4 
z e [1,3] 
z 6 [0.6,2.6] 
z £ [0.5,3] 
z~2 

local 



Beswick et al. (2008)+ 
Boyle et al. (2007)+ 
Garn fc Alexander (2009) + 



Rieke et al. (2009) 



IR; S„(24|Um) > 0.08 mjy 
IR; S v {24fim) > O.lmjy 
IR; S„(24/im) > 0.15 mjy 

IR; S„(60/«m) > 2 Jy 



Ibar et al. (2008) radio; S„(1.4 GHz) > 0.035 mjy 



Applcton et al. 

Frayer^t_aL 
^e^mour^^aL 
Yun et al. 



:2004) 

(2006) 

(2009) 

(2001) 
Garrett (2002) 
Kovacs et al. (20061 
Murphy et al. (2009b ) 
Saiina et al. (2008) 
Younger et al. (2009) 

Bell (2003) 



radio; S u (1.4 GHz) > 0.09 mjy 

radio: S„ (1.4 GHz) > 0.09 mjy 

radio: S„ (1.4 GHz) > 0.03 mjy 

radio- matched w/ S v (60 fim) > 2 Jy 
radio; 5<r WSRT sres. 
sub-mm/radio 

IR; S„(24/im) > 0.9 mjy 
IR; S v {24/im) > 0.5 mjy 

FUV (& IR) 



C"l.4GHz ~ 0.004 mjy/bm 
CTl.4 GHz ~ 0.03mJy/bm 



S„(24/xm) > 0.2 mjy 

S„(24^m) > 0.5 mjy, 
S„(70Atm) > 30 mjy 

S„(70 jum) > 15mJy 
S„(70 (im) > 6 mjy 



S„(24/mi) > 0.2 mjy 
S„ (1.4 GHz) > 0.09 mjy 



924 = 0.52-0.7 (S„(24/xm) < lmjy) 

924 = 1.39 
524 = 0.92-1.02 (S„(24|Hii) < 1 mjy) 
<?24 = 1.22 ±0.24 

(log(iTIR/£ Q ) < 11), 

g 2 4 = (-1.28 ± 0.76) + (0.22 ± 0.07) ■ log(^Tm/i ) 

(log(iTIR/L ) > 11) 

924=(0.94±0.01)-(0.01±0.01)-2 

<J24 = 0.94-l±0.25, 
970 = 2.15±0.16 

970 = 2.10 ± 0.16 (z ~ 0.2) 
570 = 2.14 ± 0.10 - 0.75 ± 0.32 • log(l + z) 
g Fm = 2.34±0.26 

<?FIR — 2 

gpiH = 2.07±0.08 

9tir = 2.41±0.39 

<3fir = 2.07/2.21 

gpiH = 2.23±0.19 

5FIR = 2.36±0.26, 
9tir = 2.64±0.26 



Radio stacking of 24 /im sources. 
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TABLE 3: 

Median of K— corrected logarithmic 24/iM/1.4GHz flux ratios, 924,0, as a function of 

REDSHIFT (CF. ALSO FlG. [131} . THE AVERAGE IR/RADIO PROPERTIES ARE GIVEN FOR THREE DIFFER- 
ENT SAMPLES OF SFGS: A RADIO-SELECTED SAMPLE, AN IR-SELECTED SAMPLE AND THE COMBINA- 
TION OF THE TWO OF THESE (FROM LEFT TO RIGHT). FOR THE JOINTLY SELECTED SAMPLE WE ALSO 
GIVE THE SCATTER IN THE RELATION (RIGHT-MOST COLUMN). ALL ERRORS ARE 2 a UNCERTAINTIES. 



redshift (924,0) (924,0) (924,0) 0-924 



radio-selected sample 


IR-selected sample 


jointly selected sample 


0.08 


< 


z 


< 





23 


1 


02« 


08 
07 


1 


35±g 


08 
12 


1 


01+0.10 
ol -0.05 


0.37±0.03 


0.23 


< 


z 


< 





33 





88i° 


16 
10 


1 


281° 


07 
07 


1 


ne+007 
ZD -0.06 


0.25±0.03 


0.33 


< 


z 


< 





45 







11 

04 


1 


291° 


04 
08 


1 


97 +0.08 
z '-0.04 


0.31±0.04 


0.45 


< 


z 


< 





67 





971° 


08 
14 


1 


261° 


05 
14 


1 


9R+0-02 
z -0.08 


0.35±0.03 


0.67 


< 


z 


< 





82 





83±° 


18 
16 


1 


291° 


06 
06 


1 


90+0. 06 
z -0.03 


0.36±0.03 


0.82 


< 


z 


< 





94 





78t° 


20 
17 


1 


321° 


03 
09 


1 


01 +0.03 
Ji -0.06 


0.35±0.03 


0.94 


< 


z 


< 


1 


14 





771° 


18 
26 


1 


361° 


07 
10 


1 


9 Q +0.07 
z -0.09 


0.49±0.04 


1.14 


< 


z 


< 


1 


51 





751°, 


12 
15 


1 


471° 


11 
07 


1 


ie+0.09 
1D -0.10 


0.91±0.06 


1.51 


< 
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TABLE 4: 

Slope (column 1) and y-Axis intercept (column 2) of the best- 
fitting LINEAR TREND TO THE EVOLUTION OF THE THREE INVESTI- 
GATED IR/RADIO PARAMETERS 924,0, 970,0 AND 9TIR- THE LINEAR 
FIT WAS PERFORMED USING THE MEDIAN LOGARITHMIC IR/RADIO RA- 
TIOS IN ALL REDSHIFT SLICES AT 2 < 1.4. THE THIRD COLUMN GIVES 
THE AVERAGE SCATTER IN THE RELATION MEASURED OVER THIS RED- 
SHIFT RANGE. THE STATES ERRORS ARE 1 a UNCERTAINTIES. 

IR/radio parameter d(q) /dz\ z <i.4. (9)2=0 (fg)U<i.4 

924 -0.015±0.136 1.275±0.098 0.417±0.038 

g 70 -0.123±0.135 2.314±0.091 0.392±0.035 

9tir -0.268±0.115 2.754±0.074 0.412±0.037 



TABLE 5: 

As for Table [3]but for the K— corrected logarithmic 70/xm/1.4GHz flux ratio 970,0- 
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TABLE 6:As for Tables [3] and [5] but for the logarithmic TIR/1.4GHz flux ratio ijtir. 
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Minimal detectable L TIR predicted based on sensitivity at: 




redshift 

Fig. 1. — : Range of total infrared luminosities Ltir (or star formation rates following iBelll (|2003[ ). see right-hand 
vertical scale) sampled between redshift and 2, given the 3cr detection limits at 24 fim, 70 (im and 1.4 GHz. The 
1.4 GHz de t ection limit has b een converted to Ltir assuming the average local TIR/radio ratio (?tir) = 2.64 (cf. 
[Belli (l200l : lYun et all (|200l - the upper and lower envelope of the arc filled with vertical light grey lines denote 
±1 a outliers to the mean of the relation). The solid tracks following a similar locus as the 1.4 GHz detection limits 
represent the smallest value of Ltir expected to be detectable based on the sensitivity at 70 /im. Different tracks 
correspond to different IR SEDs and /or libraries (s e e lege nd in lower right corner). T he thicker (solid) lines emphasize 
the predictions from the libraries of iDale fc Heloul (|2002D and iCharv fc Elbad (|200lD which are used for the IR SED 
fitting (see §[4j) . Predictions based on 3 a flux limits at 24 /iin are plotted as thin dashed tracks. An observed flux limit 
of 0.3 mJy at 24 /jm leads to a similar sampling of the IR luminosity function as is possible at 70 /im and 1.4 GHz. This 
is illustrated by the coarser dashed tracks which are joined to the fine ones with the dotted lines. Note that tracks do 
not always run across the entire redshift range due to the finite choice of IR luminosities represented in the various 
SED libaries. 
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radio— selected sample: 



2- band matching w/ 34 fim 
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d(1.4 GHz, counterpart) [arcsec/max(FWHM 1 4 GHz ,FWHM band 2 )] 




24 /xm-selected sample: 

band matching w/ 1.4 GHz 
band matching w/ 70 /xm 
band matching w/ optical 

" " " band matching w/ IRAC 



(b) 



0.1 0.15 0.3 0.35 0.3 0.35 0.4 0.45 0.5 

d(24 fj,m, counterpart) [arcsec/max(FWHM 24 ^ m >FWHM band 2 )] 

Fig. 2. — : (a) - Radio-selected sample; separation between radio source position and nearest catalog counterpart at 
24 (black line), 70 /xm (light grey line) and at optical wavelengths (dark grey). The distance is expressed in fractions 
of the FWHM of the larger PSF of the two bands involved in the match. (The resolution is 0.6", 2.5", 5.8" and 18.6" 
in the optical and at 1.4 GHz, 24 /im and 70 /im, respectively. The IRAC PSF has a FWHM which ranges from 1.6" to 
2" between the 3.6 /xm and the 8 /im channel. For the histogram in the lower panel of the figure we assume an average 
FWHM of 1.8".) 

(b) - IR-selected sample; separation (defined as in (a)) between 24 ^m source position and the corresponding counterpart 
at 1.4 GHz (black line), 70 ^m (light grey line) and at optical (solid dark grey line) and NIR (IRAC) wavelengths (dotted 
dark grey histogram). 




<S/N 70 > 



10 



50 



100 



1000 
500 E- 



w 100 r 
| 50 §■ 

o 
o 

10 t- 

5 L 



-i 1 r 



-i 1 r- 



b) 



nr 



. 1 1 1 1 I 1 

■ PRF flux 

□ upper lim. from local noise 




10 50 
S„(70Mm) [mJy] 



100 



Fig. 3.—: (a) - Histogram showing the available flux information for the 24 /xm counterparts to radio-selected COSMOS 
sources in a given bin of 24 /zm flux density (grey area - PSF-flux measurement with S/N > 3; white area - 3 a upper 
limit from local noise). 

(b) - Same information as displayed in upper panel but for flux constraints on 70 /im counterparts of radio-selected 
sources (grey area - PRF-flux measurement with S/N > 3; white area -3a upper limit from local noise). 
The 24 /im flux distribution shows a very sharp transition from PSF-fitted detections to upper limits. This is due to 
the fact that the 24 ^m catalog is a flux- limited catalog while the 70 /im catalog is selected according to S/N. 
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redshill redshlft 

Fig. 4. — : Bar diagrams showing the quality of the available flux information at different rcdshifts for sources in the 
radio-selected (left) and IR-sclcctcd sample (right). 

Left - Dependence on rcdshift of the fraction of radio-selected sources (referred to the total number of objects in a 
given redshift slice) with a positive identification in both the 24 /xm and 70 /xm catalog (black histogram) , as well as 
the fractional contribution of sources with flux limits in cither one of the two bands (light grey - detected 24 /tm 
counterpart, upper flux limit at 70 /xm; dark grey - upper flux limit at 24 /xm, detection at 70 /xm) or in both (white 
area) . 

Right - As for the first panel but for the IR-selected sample (restricted to sources with S'„(24 /xm) > 0.3 mJy). Black 
- counterpart detected at both 70 /xm and 1.4 GHz; dark grey - detected counterpart at 70 /xm, upper flux limit at 
1.4 GHz; light grey - detected counterpart at 1.4 GHz, upper flux limit at 70 /xm; white - upper flux limits at both 
70 /xm and 1.4 GHz. 
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Fig. 5.—: (a) - Histogram displaying the contribution of well-defined flux measurements and upper flux limits for the 
70 /im counterparts of S-COSMOS sources (with S U (2A /im) > 0.3 mJy) selected at 24 /xm as a function of their 70 /im 
flux density (grey area - PRF-flux measurement with S/N > 3; white area - 3 a upper limit local noise). 
(b) - Same information as in the upper panel of the figure but for flux constraints at 1.4 GHz (grey area - flux 
measurement with S/N > 3; white area - 3 a upper limit from local noise). 
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Fig. 6. — : Location of sources selected either only at 24 /.tm (orange symbols) or 1.4 GHz (green) or present in both 
catalogs (black) in plots comparing fluxes that will be used for the computation of their IR/radio properties: (a) - 
70 fim flux constraints S v (70 fj,m) as a function of 1.4 GHz flux density Sj, (1.4 GHz); (b) - S v (70 fim) vs. 5„(24 fim); 
(c) - S„(24fim) vs. ^(1.4 GHz). 

The IR-selected sample comprises sources plotted in orange and black. The union of green and black symbols defines 
the radio-selected sample (see § 12. 3|) . In panel c) the straight cut-offs confining the distribution of the green points on 
the left and that of the orange points on the bottom reflect the lower flux limits of the radio- and IR-selected sample, 
respectively. 
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(u-K) 

Fig. 7. — : Correlation between the rest frame colours PI and (u—K) for radio-selected sources from the VLA-COSMOS 
survey. The bisector fit - see equation ([3]) - to the relation is shown with a b l ack lin e; grey lines illustrate the uncertainty 
in the fit. The separation between SF and AGN galaxies of Smolcic ct al. (2008) at Pl= 0.15 corresponds to a cut at 
(u — K) — 2.42 (cf. dotted line). Note that due to our treatment of composite SF/AGN sources we adopt a slightly 
different colour threshold for the selection of SF systems (see following figure). Upper left corner: histogram of offsets 
API (measured parallel to the vertical axis for each data point) from the black trend line. The dispersion of the 
relation is ctapi = 0.17. 
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(u-K) 

Fig. 8. — : (a) - Distr ibution with respect to PI colour of star forming (SF), AGN and composite systems in the 
SDSS/NVSS sample of lSmolcic etlll (f2008h . (b) - Probability Pr cff . (SF | PI) (cf. equation ©) of correctly classifying 
a galaxy as star forming (black line) or AGN (dark grey line) at a given value of PI colour based on panel (a) of this 
figure, (c) — Probability Pr e ff. (SF) (black curve), expressed as a function of rest frame (u — K) colour. The conversion 
from panel (b ) is based on the average trend and the scatter of the relation shown in Fig. [Jj The adjacent light grey area 
illustrates the error on the Pr c ff. (SF | (u — K)) curve induced by the uncertainty of the best-fit linear relation between 
PI and (it — K). In dark grey the complementary probability of correct AGN-classification, Pr e ff. (AGN | (u — K)), 
is plotted (without the error which mirrors the one of the curve for SF systems). Equal probabilities of correctly 
classifying objects as SF or AGN, respectively, are reached at (u — K) = 2.36. 
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Fig. 9. — : Distribution of rest frame (u — K) colours of galaxies in the radio-selected sample (light grey) and the 
IR-selected sample (limited to sources with S' 1/ (24/xto) > 0.3 mJy; dark grey) in six different redshift slices. The 
hatched areas indicate those sources which are common to both samples. The median (u — K) colours in the radio- 
and IR-selected samples are marked by the vertical light and dark grey lines, respectively. Along the upper edge of 
the plot the probability Pr (SF | (u ~ K j) of 'SF-hood' associated with a given value of (u — K) colour is indicated. 
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Fig. 10. — : Schematic illustration of the origin of the difference Agbi as between the average IR/radio ratio measured 
using an IR- or a radio-selected sample (see ij l5.ll for details). The vertically offset groups of curves show typical 
IR-radio SEDs of three objects with different observed bolometric flux. The increasing intensity (top to bottom) of the 
colour grey with which each of the three groups is drawn reflects that source counts rise with decreasing flux according 
to dN /ds oc S~ . Within each group of curves with a similar observed flux density, the central one represents the SED 
of an object with an average IR/radio ratio. The dashed (dotted) curve indicate +3 a (-3 a) outliers to the IR-radio 
relation. Arbitrary observational limits in the radio (left) and MIR (right) window are marked with horizontal black 
bars. 
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Fig. 11. — : Ratio of observed (i.e. not ^-corrected) 24 fim and 1.4 GHz fluxes, ^24, for galaxies classified as SF (top) 
and AGN (bottom). Sources with the highest probability of being star forming are plotted in dark blue, while those 
with the highest probability of being AGN are shown in red (see lower left corner of both panels). The dotted vertical 
lines mark the redshifts above which only (U)LIRGs remain in the sample. Large dots are used for sources that are in 
both the IR- and radio-selected sample; objects plotted with smaller dots are found in only one of the two samples. 
The style of the symbol reveals if a source has a spectroscopically (filled dot) or photometrically (open circle) measured 
redshift. This distinction between the quality of the redshift information is not made for lower and and upper limits 
(from the IR- and radio-selected sample, respectively) on ^24 which are reported with fainter colours. 
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Fig. 12. — : As for Fig. [IT] but for the ratio of observed 70 /im and 1.4 GHz flux, (770. obs- Symbols and colours are 
identical to those used in Fig. QTJ 
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Fig. 13. — : Evolution of the if-corrected 24/zm/1.4GHz flux ratios, 924, 0, as a function of redshift. Upper panel - 
Cumulative/probability distribution functions of 924,0 in a radio- (light grey) and an IR-selected (dark grey) sample 
of SF galaxies, as well as in the union of the two (black curves). The panels shown here represent every second bin of 
a finer redshift sampling presented in full in the lower panel of the figure. The curves give the probability of finding a 
system with a smaller (or - in the case of the IR-selected sample - a larger) 924,0 than a given value on the horizontal 
axis. Hatched areas span the 95% confidence interval of the radio- and IR-selected distribution functions at fixed 924,0- 
In each redshift bin the dashed red line traces the Gaussian distribution which fits the measured distribution function 
of the joint IR- and radio-selected samples best. The intersection of the black curve with the 50% probability line 
(dotted horizontal line) defines the median value of 924,0 in the three different samples (see also Table [3]). 
Lower panel - if-corrected 24/im/1.4GHz flux ratios, 924,01 as a function of redshift. Symbol colours and style are 
identical to those used in Fig. [TTJ The black dots and error bars mark the medians and associated 2 a errors obtained 
from the distribution function of 924, of the joint IR- and radio-selected sample (see black curves in the upper panel). 
Width and location of the redshift bins were chosen such as to always include the same number of objects. The black 
line represents the best-fitting evolutionary trend (fit as a function of linear redshift space) to the medians at z < 1.4. 
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Fig. 14. — : Variation as a function of the threshold used to select SF systems of: (top) - the extrapolation to z = 
of the evolutionary trend line fitted to the medians (924,0} at z < 1.4; (centre) - the mean scatter in the population, 
averaged over all redshift slices at z < 1.4; (bottom) - the slope of the evolutionary trend of (724,0 at z < 1.4. 
In each panel the black point highlights the measurement obtained if all objects with probability Pr (SF) > 50% are 
considered star forming (the convention used when plotting Figs. fTTI and IT3| . The scatter decreases with increasing 
purity of the star forming sample, however no significant changes in the best-fitting evolutionary trend line are seen 
when samples consisting of sources with a probability of at least 20%, 30%, ... 80% of being SF are constructed. 
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Fig. 15. — : As for Fig. [13] but showing the distribution functions of the -ftT-corrected IR/radio flux ratio 570,0 (top) as 
well as the redshift evolution of the median ratio (970,0) at z < 1.4 (bottom). 
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Fig. 16. — : As for Fig. [14] but investigating variations in the scatter and evolutionary trend of 970.0 if different 
thresholds for the selection of a sample of SF galaxies are adopted. 
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Fig. 17. — : As for Figs. [T3l and fT5l but showing the distribution functions of the TIR/radio flux ratio qtir (top) as well 
as the redshift evolution of the median ra tio (<7tir) at z < 1-4 (botto m). The vertical dashed line marks the locally 
measured average TIR/radio ratio of 2.64 (jYun et all [2001; Belli l200l . 
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Fig. 18. — : As for Figs. [T4l and [T6l but assessing the robustness of the best-fitting trend for the redshift evolution of 
9tir to changing the probability threshold for the selection of a sample of SF galaxies. Also shown are changes in the 
scatter of the TIR/radio correlation (middle) as a function of varying sample selection. 
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Fig. 19. — : Comparison of the IR/radio properties of model SF systems (coloured tracks; see legend along upper 
edge) with those of AGN-bearing sources detected in the XMM-Newton survey of the COSMOS field, (a) - observed 
24/mi/1.4GHz flux density ratio 924, obs! (b) - observed 70/mi/1.4GHz flux density ratio 970, obs- 
Large symbols are used for sources that are found in both the IR- and radio-selected sample, small dots if a given 
source is only found in one of the two samples. 
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Fig. 20. — : Assessment of the contribution of AGN-bearing sources, fAGN (924, obs) , to the total number of objects lying 
on the IR/radio locus expected for SF systems. 

Upper panel — Definition of the redshift-dependent locus of SFGs (delimited by the solid black lines marked with 
'924, obs — (924, template) = ±0.7', respectively; see text for details) in which the relative frequency of AGN and SF 
sources is mapped. The coloured tracks show the evolution of the observed 24/zm/1.4GHz flux ratios, 924. obs, for 
different IR-SEDs from a variety of template libraries. The vertical dotted lines indicate the limits of the redshift bins 
used in the lower part of this figure. 

Lower panel - Variation of the relative contribution of AGN and SF sources, fAGN, to the total population at a given 
value of (?24,obs- The changes are traced between ±2er of the expected mean - (924, template) U _ for local IR SEDs 
at different redshifts (see upper panel of the figure and details in text). The black line is a smoothed version of 
fAGN (924, obs) , obtained by taking a 5-point running average of the finer and noisier mapping reported in light grey. 
The hatched area shows the associated ±2cr uncertainty region. 
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Fig. 21. — : As for Fig. [20]but measuring the relative frequency fAGN(<?7o,obs) of AGN and SF galaxies in a ±2cr band 
(cf. black lines in upper panel) around the mean expected value of the observed 70 ^m/1.4 GHz flux ratio (gy p, te mplate) \z 
based on template SEDs of star forming galaxies. Symbols and colours are identical to those used in Fig. [^01 
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Fig. 22. — : Dependence of the TIR/radio flux ratio (?tir of star forming galaxies on total IR luminosity I/tir (top) 
and 1.4 GHz luminosity £i.4gh~z (bottom). Medians and their associated errors derived with survival analysis are 
plotted in black. In the thinner two windows the scatter in the measurements of giro is shown for different luminosity 
bins. A double off-set binning scheme was chosen in order to account for the uncertain x-axis value of upper limits in 
luminosity (diagonal arrows). 
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Fig. 23. — : TIR/radio ratios of all IR- or radio-selected sources in the sample that lie at z > 2.5 (sources are coloured 
according to their probability of 'SF-hood'; see colour scale). On the right hand side of the plot the cumulative 
distribution function is plotted together with the best- fit Gaussian distribution (cf. Fig. IT?]) . The dashed line marks 
the locally measured average TIR/radio ratio. 



